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Research  Center.  Dr.  Robert  W.  Whalin,  Technical  Director,  WES,  was  the 
director  of  the  study;  and  Mr.  D.  L.  Robey,  Chief  of  the  Ecosystem  Research 
and  Simulation  Division  (ERSD),  Environmental  Laboratory,  was  the  Study 
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ED;  and  Mr.  Butler.  Dr.  Y.  P.  Sheng  of  University  of  Florida  was  a  consultant 
in  the  early  phases  of  the  study. 

The  June-July  1980  data  presented  in  the  report  were  obtained  from 
Dr.  W.  Boicourt  of  the  Horn  Point  Environmental  Laboratories  in  Cambridge,  MD. 
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This  report  was  edited  by  Mrs.  M.  C.  Gay,  Information  Technology 
Laboratory,  WES. 

Commander  and  Director  of  WES  during  preparation  of  this  report  was 
COL  Larry  B.  Fulton,  EN.  Technical  Director  was  Dr.  Robert  W.  Whalin. 


The  report  should  be  cited  as  follows: 


Joimson,  B.  H.,  Kim,  K.  W. ,  Heath,  R.  E. ,  Hsieh, 
H.  L.  1991.  "Development  and  Verification  of  a 
Numerical  Hydrodynamic,  Salinity  and  Temperature 
Bay,"  Technical  Report  HL-91-7,  US  Army  Engineer 
Station,  Vicksburg,  MS. 


B .  B . ,  and  Butler , 
Three-Dimens ional 
Model  of  Chesapeake 
Waterways  Experiment 


1 


CONTENTS 


Page 

PREFACE .  1 

CONVERSION  FACTORS,  NON-SI  TO  SI  (METRIC) 

UNITS  OF  MEASUREMENT .  4 

PART  I:  INTRODUCTION .  5 

Background .  5 

Purpose .  6 

Scope .  6 

PART  II:  MODEL  DESCRIPTION .  9 

Boundary-Fitted  Equations .  9 

Boundary  Conditions .  12 

Numerical  Solution  Algorithm .  14 

Turbulence  Parameterization .  15 

PART  III;  APPLICATION  WITH  AN  EARLY  VERSION  OF  THE  MODEL  TO  THE 

JUNE-JULY  1980  AND  APRIL  1983  DATA  SETS .  20 

Strategy .  20 

June-July  1980 .  21 

April  1983 .  24 

Cone  lus  i  ons .  26 

PART  IV:  VERIFICATION  OF  THE  FINAL  MODEL  TO  THE 

SEPTEMBER  1983  DATA  SET .  27 

Description  of  Field  Data .  27 

Initial  Conditions .  27 

Boundary  Forcing  Data .  28 

Model  Coefficients .  28 

Verification  Results .  29 

Conclusions .  32 

PART  V:  MODEL  VERIFICATION  ON  AN  ANNUAL  TIME  SCALE .  33 

1984  Simulation .  33 

1985  Simulation .  36 

1986  Simulation .  38 

PART  VI:  SUMMARY  AND  CONCLUSIONS .  40 

Summary .  40 

Conclusions .  41 

REFERENCES .  43 

TABLES  1-15 

FIGURES  1-164 

APPENDIX  A*:  1984  RESULTS .  Al 


*  A  limited  number  of  copies  of  Appendixes  A  through  C  are  published  under 
separate  cover.  Copies  are  available  from  National  Technical  Information 
Service,  Springfield,  VA  22161. 


2 


Page 


APPENDIX  B*:  1985  RESULTS .  B1 

APPENDIX  C*;  1986  RESULTS .  Cl 

APPENDIX  D:  NOTATION .  D1 


A-l 


3 


CONVERSION  FACTORS,  NON-SI  TO  SI  (METRIC) 
UNITS  OF  MEASUREMENT 


Non— SI  units  of  measurement  used  in  this  report  can  be  converted  to  SI 
(metric)  units  as  follows: 

_ Multiply _  _ _  To  Obtain 

cubic  feet  0.02831685  cubic  metres 

feet  0.3048  metres 


DEVELOPMENT  AND  VERIFICATION  OF  A  THREE-DIMENSIONAL  NUMERICAL 


HYDRODYNAMIC,  SALINITY,  AND  TEMPERATURE  MODEL  OF  CHESAPEAKE  BAY 

PART  I :  INTRODUCTION 

Background 

1.  The  Chesapeake  Bay  is  one  of  the  Nation's  most  valuable  natural 
resources.  It  supports  important  commercial  and  recreational  fisheries, 
transportation,  industry,  recreation,  and  tourism,  and  provides  irreplaceable 
habitat  for  marine  resources  and  wildlife.  However,  the  estuary  has  been 
subjected  to  increasing  environmental  stress  in  recent  decades,  and  the  pro¬ 
ductivity  and  beauty  of  the  Chesapeake  Bay  have  significantly  declined.  In 
1983,  the  US  Environmental  Protection  Agency  (USEPA)  identified  major  contri¬ 
buting  factors  to  the  bay's  decline  as  inputs  of  nutrients  and  toxicants  from 
point  and  nonpoint  sources,  changes  in  land  use  within  the  basin  with  result¬ 
ing  modification  of  the  watershed  and  resource  habitat,  and  concurrent  impacts 
of  natural  events  such  as  floods  and  droughts  (USEPA  1983a,  1983b).  Because 
population  within  the  bay  drainage  basin  is  still  increasing,  and  development 
pressures  will  continue  into  the  foreseeable  future,  it  is  necessary  that 
strategies  are  developed  to  reverse  the  present  bay  decline  and  to  accommodate 
future  growth  in  an  environmentally  sound  manner. 

2.  The  Chesapeake  Bay  Program  (CBP) ,  established  in  1983,  provides  a 
management  structure  through  which  the  activities  of  state  and  Federal  agen¬ 
cies,  as  well  as  those  of  private  citizens,  can  be  coordinated  toward  the  goal 
of  bay  restoration.  A  number  of  tools  have  been  developed  by  the  CBP  that 
will  assist  in  the  planning,  implementation,  and  evaluation  of  these  strate¬ 
gies.  Among  these  tools  are  the  following;  a  baywide  long-term  monitoring 
program  to  determine  water  quality  conditions  and  trends;  a  comprehensive 
Chesapeake  Bay  data  base;  and  a  series  of  numerical  models  to  evaluate  alter¬ 
native  control  strategies  and  to  guide  the  establishment  of  pollutant  reduc¬ 
tion  goals.  In  June  1985,  the  CBP  Implementation  Committee  approved  a 
modeling  strategy  that  called  for  phased  development  of  these  models.  Phase  I 
was  the  refinement,  computer  code  conversion,  and  updating  of  the  existing 
Chesapeake  Bay  Watershed  Model,  which  predicts  the  delivery  of  nutrients  to 
the  estuary  from  point  and  nonpoint  sources  above  the  fall  line  and  nonpoint 
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sources  below  the  fall  line.  Phase  II  was  the  development  of  a  s teady-s te t c 
(coarse  grid)  water  quality  model  of  the  bay  to  as.sess  general  response  of  trie 
system  to  nutrients  and  the  relative  importance  of  various  processes.  The 
third  phase  of  the  modeling  strategy  consisted  of  a  three-dimensional  (3-D;, 
time-varying,  hydrodynamic,  and  water  quality  model  for  the  bay  and  tribu¬ 
taries  that  could  provide  a  detailed  assessment  of  the  system's  response  to 
nutrient  inputs  and  other  parameters  varied  realistically  over  time  and  space. 
The  3— D  model  represents  the  means  through  which  proposed  management  actions 
call  be  tested  before  implementation,  allowing  more  cost-effective  selection  of 
appropriate  strategies. 


Purpose 

3.  It  has  been  acknowledged  that  water  quality  impacts  in  the  Chesa¬ 
peake  Bay  cannot  be  successfully  assessed  without  an  accurate  description  of 
the  hydrodynamic  processes.  Thus,  the  purpose  of  this  study  was  the  develop¬ 
ment  and  subsequent  vei if icar ion  of  a  3-D  hydrodynamic  model  of  Chesapeake 
Bay.  The  particular  hydrodynamic  computer  code  employed  is  called  CH3D  (Cur¬ 
vilinear  Hydrodynamics  in  Three  Dimensions) .  The  basic  code  was  developed  by 
Sheng  (1986)  for  the  US  Army  Engineer  Waterways  Experiment  Station  (WES)  but 
has  been  extensively  modified  in  its  application  to  Chesapeake  Bay.  As  its 
name  implies,  CH3D  makes  hydrodynamic  computations  on  a  curvilinear  or 
boundary-fitted  planform  grid.  Physical  processes  impacting  baywide  circula¬ 
tion  and  vertical  mixing  that  are  modeled  include  tides,  wind,  density  effects 
(salinity  and  temperature),  freshwater  inflows,  turbulence,  and  the  effect  of 
tha  earth's  rotation.  Adequately  representing  the  vertical  turbulence  is 
crucial  to  a  successful  simulation  of  stratification,  destratification,  and 
anoxia  in  the  bay.  The  boundary-fitted  coordinates  feature  of  the  model  pro¬ 
vides  enhancement  to  fit  the  deep  navigation  channel  and  irregular  shoreline 
configuration  of  the  bay  and  permits  adoption  of  an  accurate  and  economical 
grid  schematization. 


Scope 

4.  The  Chesapeake  Bay  is  one  of  the  largest  estuaries  in  the  world.  As 
shown  in  Figure  1,  the  main  bay  extends  approximately  300  km  north  from  the 
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ocean  entrance  to  the  Susquehanna  River.  The  average  depth  of  the  bav  is 
about  8  m,  although  a  natural  channel  with  depths  greati-r  than  I'd  m  tia\'ersc.s 
the  bay  for  more  than  60  percent  of  its  length,  Tiie  bay  i.s  iiregular  in 
shape,  varying  in  width  from  6.4  km  between  Annapolis,  .MD ,  and  Kent  Island,  in. 
the  Upper  Chesapeake  Bay,  to  48.3  km  in  the  middle  bav  off  the  Potomac  River. 
The  bay  is  long  enough  to  accommodate  one  complete  tidal  wave  at  all  times. 

5.  The  numerical  grid  employed  in  the  3-D  hydrodynamic  model  is  shown 
in  Figure  2.  There  are  734  active  horizontal  cells  and  a  maximum  of  15  verti¬ 
cal  layers,  resulting  in  over  4,000  computational  cells.  To  capture  the 
important  features  of  the  hydrodynamic  processes  and  bathymetry  in  the  bay, 
grid  resolution  is  1.52  m  vertical  and  approximately  10  km  longitudinal  and 

3  km  lateral.  Major  tributaries,  i.e.,  the  James,  York,  Rappahannock, 

Potomac,  Patuxent  and  Susquehanna  Rivers,  are  modeled  fully  3-D  in  the  lower 
reach  and  two-dimensional  (2-D)  in  the  upper  reach  with  a  single  cell  spanning 
the  width  of  the  tributary  Based  upon  an  average  inflow  of  70,000  cfs,  typi¬ 
cally  about  90  percent  of  the  fresh  water  in  the  bay  enters  by  these  rivers. 

6.  A  successful  verification  of  the  hydrodynamic  model  requires  sets  of 
synoptic  data.  These  sets  must  contain  freshwater  inflows  and  their  tempera¬ 
ture  on  the  major  tributaries;  tides  at  the  bay  entrance  as  well  as  at  various 
interior  stations;  meteorological  data  at  one  or  more  stations  from  which  the 
surface  wind  stress  and  heat  flux  can  be  determined;  and  currents,  tempera¬ 
ture,  and  salinity  at  several  locations  throughout  the  bay  as  well  as  at  the 
ocean  boundary. 

7.  The  bay  being  so  large  contributes  to  the  fact  that  there  is  a  lack 
of  synoptic  data  throughout  the  bay  and  its  tributaries.  Three  relatively 
extensive  synoptic  data  sets  were  identified  for  .se  in  the  initial  verifica¬ 
tion  of  the  3-D  model.  These  data  sets  were  collected  during  June-July  1980, 
April  1983,  and  September  1983.  The  1980  data  were  collected  and  provided  by 
Dr.  William  Boicourt  of  the  Horn  Point  Environmental  Laboratories  at  Cam¬ 
bridge,  MD,  and  the  1983  data  were  collected  primarily  by  the  National  Ocean 
Service  (NOS)  and  provided  by  Dr.  Alan  Blumberg  of  HydroQual ,  Inc,  Mahwah,  NJ . 
The  1980  data  set  represents  low  inflow  conditions  and  is  a  good  characteriza¬ 
tion  of  the  summer  circulation  that  may  occur  in  any  year.  The  spring  1983 
data  represent  conditions  during  a  large  spring  runoff  event,  whereas  the  fall 
1983  data  set  contains  a  strong  wind-mixing  event  that  resulted  in  a  destrati¬ 
fication  of  the  bay.  This  process  was  aided  by  a  rapid  cooling  of  the  surface 
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waters.  Final  verification  of  the  3-D  model  was  achieved  through  yearlong 
simulations  of  the  years  1984,  1985,  and  1986,  These  years  represent  high, 
low,  and  average  freshwater  inflows,  respectively.  Table  1  summarizes  condi¬ 
tions  existing  during  each  of  these  simulation  periods, 

8.  Results  from  each  of  these  six  simulations  are  presented.  It  should 
be  noted  that  results  from  the  June-July  1980  and  April  1983  simulations  are 
from  the  application  of  an  earlier  version  of  the  final  model.  Thus,  these 
results  should  be  viewed  only  in  a  qualitative  s-^nse.  Application  to  the  fall 
1983  data  is  discussed  in  detail  since  this  data  set  guided  the  modifications 
resulting  in  the  model  that  was  ^nally  applied  to  the  yearlong  production 
simulations.  These  yearlong  simulations  also  served  as  additional  verifica¬ 
tion  of  the  3— D  model  since  no  additional  adjustment  of  model  coefficients  was 
made.  Only  selected  results  from  the  yearlong  simulations  are  presented. 
Complete  results  are  given  in  Appendixes  .\-C. 

9.  Before  discussing  verification  results,  a  limited  discussion  of  the 
numerical  model  is  provided.  For  more  information,  see  Johnson  et  al .  (1991). 
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PART  II:  MULEL  DESCRY  TION 


10.  To  better  resolve  complex  geonietrie-;  in  the  horizontal  directions, 
cr.c  Che  .apeake  Bay  3-D  model  makes  computations  on  the  boundary-fitted  or  gen¬ 
eralized  curvilinear  planforra  grid  shown  in  Figure  2.  This  necessitates  the 
transformation  of  the  governing  equations  into  boundary-fitted  coordinates 

If  only  the  (x,y)  coordinates  are  transforined ,  a  system  of  equations 
similar  to  those  solved  by  Johnson  (1980)  for  vertically  av  raged  i iow  fields 
is  obtained.  However,  in  the  Chesapeake  Bay  model,  not  only  are  the  (x,y) 
coordinates  transformed  into  the  ,r)'  curvilinear  system  but  also  the 
velocity  is  transforme.,  su^n  that  its  components  are  perpendicular  to  the 
(^.^)  coordinate  lines.  This  is  accomplished  by  employing  the  following  defi¬ 
nitions  for  the  components  of  the  Cartesian  'elocity  (u,v)  in  terms  of 
cont ravar iant  components  u  and  v 


u  =  x^u  +  x^v 


(1) 


V  =  y^u  I-  y,v  (2) 

With  tl  9  governing  equatio.is  written  in  terms  of  the  contravariant  components 
of  the  ve’-icity,  boundary  conditions  can  be  prescribed  on  a  boundarv-f itted 
grid  in  the  same  manner  as  on  a  Cartesian  grid  since  u  and  v  are  perpen¬ 
dicular  to  the  curvilinear  cell  faces;  e.g.,  at  a  land  boundary,  either  u 
or  v  is  set  to  zero. 

11.  Initially  the  vertical  dimension  was  handled  through  the  use  of 
what  is  commonly  called  a  sigma  stretched  grid.  However,  with  the  grid  reso¬ 
lution  shown  in  Figure  2,  it  was  observed  that  stratification  in  the  deep 
channels  could  not  be  maintained  during  long-term  simulations  on  Chesapeake 
Bay.  With  a  sigma  stretched  grid,  the  bottom  layer  in  one  column  communicates 


*  For  c  nvenience,  symbols  and  unusual  abbreviations  are  listed  and  defined 
in  the  Notation  (Appendix  D) . 
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with  the  bottom  layer  in  an  adjacent  column.  Thus,  if  depth  changes  are 
rather  coarsely  resolved,  channel  stratification  cannot  be  maintained.  As  a 
result,  the  governing  equations  were  rederived  for  solution  on  the  Cartesian 
or  z— plane  in  the  vertical  direction. 

12.  With  both  the  Cartesian  coordinates  and  the  Cartesian  velocity 
transformed,  the  following  boundary— fitted  equations  for  u  ,  v  ,  w  ,  S  , 
and  T  to  be  solved  in  each  vertical  layer  are  obtained: 
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where 

G22  =  ^  yv 

Gi2  =  Gzi  =  +  y^y^ 

*^11  ~  *  yf 

Equations  3-7  have  been  written  in  terms  of  dimensionless  variables  where 
h  =  layer  thickness 
t  =  time 

5"  -  water-surface  elevation 
w  =  vertical  component  of  velocity 
p  -  water  density 
S  -  salinity 
X  =  temperature 

Xhe  dimensionless  parameters  Ro  .  E^  •  Etp  ,  Fr  ,  and  Pr^  are  defined  as 

Ro  =  u,/fx. 
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Ev  = 

Fr^,  =  Fi-/[(^,  - 
Fr  =  Uy(gZ,)i/2 

Prv  = 


where 

A^,K^  “  vertical  turbulent  eddy  coefficients 
f  —  Coriolis  parameter 
g  —  acceleration  due  to  gravity 

and  the  quantities  ,  and  are  arbitrary 

reference  values  of  the  velocity,  density,  horizontal  and  vertical  dimensions, 
vertical  viscosity,  and  vertical  diffusivity,  respectively.  An  equation  of 
state  relating  the  water  density  to  the  salinity  and  temperature  closes  the 
system. 

13.  As  noted,  there  are  a  maximum  of  15  vertical  layers  in  the 
Chesapeake  Bay  model.  Each  is  1.52  m  thick  except  for  the  top  layer,  which 
varies  with  the  tide.  The  water— surface  slope  terms  in  Equations  3  and  4  are 
evaluated  from  similar  transformed  equations  for  the  vertically  averaged  flow 
field.  The  horizontal  diffusion  terras  are  quite  lengthy  and  thus  are  not 
listed  here,  but  are  included  in  Johnson  et  al .  (1991).  It  might  be  noted 
that  initially  the  convective  terms  in  the  momentum  equations  resulted  in 
unstable  solutions  due  to  being  written  with  velocity  squared  terms  present 
(see  Sheng  1986).  However,  writing  these  terms  as  presented  in  Equations  3 
and  4  solved  this  problem. 


Boundary  Conditions 


14.  The  boundary  conditions  at  the  free  surface  are 


r 


(8) 


12 


(9) 


ax 


"  K  (T  -  TJ 


4^  ==0 


(10) 


whereas  the  boundary  conditions  at  the  bottom  are 
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where 

-  components  of  surface  shear  stress 
C  -  surface  drag  coefficient 

W  -  wind  speed 

K  -  surface  heat  exchange  coefficient 
X,  -  equilibrium  temperature 

-  components  of  bottom  shear  stress 
Cj  -  bottom  friction  coefficient 

Ui,Vi  -  horizontal  velocity  components  next  to  the  bottom 
With  Zj  equal  to  one-half  the  bottom  layer  thickness,  is  given  by 


C<i  =  k2  {[In  (Zi/z„)]) 


(13) 


where  k  is  the  von  Karman  constant  and  z.„  is  the  bottom  roughness  height. 
As  presented  by  Garratt  (1977) ,  the  surface  drag  coefficient  is  computed  from 


C  =  (0. 75  +  0.067  W)  X  lO'^ 


(1^) 
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with  the  maximum  allowable  value  being  0.003.  The  surfai’e  heat  exchange  coef¬ 
ficient  K  and  the  equilibrium  temperature  T,  are  computed  from  meteorolog¬ 
ical  data  as  discussed  by  Edinger,  Brady,  and  Geyer  (1974). 

15.  Along  the  shoreline  where  river  inflow  occurs,  the  freshwater  in¬ 
flow  and  its  temperature  are  prescribed  and  the  salinity  is  assumed  to  be 
zero.  At  an  ocean  boundary,  the  water— surface  elevation  is  prescribed  along 
with  time-varying  vertical  distributions  of  salinity  and  temperature.  During 
flood,  the  specified  values  of  salinity  and  temperature  are  employed,  whereas 
during  ebb,  interior  values  are  advected  out  of  the  grid.  Along  a  solid 
boundary,  the  normal  component  of  the  velocity  and  the  viscosity  and 
diffusivity  are  set  to  zero. 

Numerical  Solution  Algorithm 

16.  Finite  differences  are  used  to  replace  derivatives  in  the  governing 
equations,  resulting  in  a  system  of  linear  algebraic  equations  to  be  solved. 

As  discussed  by  Johnson  et  al .  (1991),  both  external  and  internal  mode 
equations  are  solved. 

17.  The  external  mode  consists  of  equations  for  the  water-surface  ele¬ 
vation  C  and  vertically  integrated  contravariant  unit  flows  U  and  V  . 

All  of  the  terms  in  the  transformed  vertically  averaged  continuity  equation 
are  treated  implicitly,  whereas  only  the  water-surface  slope  terms  in  the 
transformed  vertically  averaged  momentum  equations  are  treated  implicitly. 
Those  terms  treated  implicitly  are  weighted  between  the  new  and  old  time- 
steps.  The  resulting,  finite  difference  equations  are  then  factored  such  that 
a  ^-sweep  followed  by  an  fj-sweep  of  the  horizontal  grid  yields  the  solution  at 
the  new  time— step. 

18.  The  internal  mode  consists  of  computations  from  Equations  3—7  for 
the  three  velocity  components  u  ,  v  ,  and  w  ;  salinity;  and  temperature. 

The  only  terms  treated  implicitly  are  the  vertical  diffusion  terms  in  all 
equations  and  the  bottom  friction  and  surface  slope  terms  in  the  momentum 
equations.  Values  of  the  water-surface  elevations  from  the  external  mode  are 
used  to  evaluate  the  surface  slope  terms  in  Equations  3  and  4.  As  a  result, 
the  extremely  restrictive  speed  of  a  free-surface  gravity  wave  is  removed  from 
the  stability  criteria.  Upwind  differencing  is  used  to  represent  the  convec¬ 
tive  terms  in  the  momentum  equations,  whereas  a  spatially  third-order  scheme 


14 


developed  by  Leonard  (1979)  called  QUICKEST  is  used  to  represent  the  advective 
terms  in  Equations  6  and  7  for  salinity  and  temperature. 

19.  It  should  be  noted  that  once  the  u  and  v  velocity  components 
are  computed,  they  are  slightly  adjusted  to  ensure  the  conservation  of  mass. 
This  is  accomplished  by  forcing  the  sum  of  u  over  the  vertical  to  be  the 
vertically  averaged  velocity  U/H  and  the  sum  of  v  over  the  vertical  to 
equal  V/H  .  where  H  is  the  total  water  depth. 

20.  A  staggered  grid  is  used  in  both  the  horizontal  and  vertical  direc¬ 
tions  of  the  computational  domain.  In  the  horizontal  directions,  a  unit  cell 
consists  of  a  ("-point  in  the  center,  a  u-point  to  its  left,  and  a  v-point  to 
its  bottom.  In  the  vertical  direction,  the  vertical  velocities  are  computed 
on  the  upper  and  lower  cell  faces.  Temperature,  salinity,  and  density  are 
computed  at  the  center  of  the  3-D  cell. 

Turbulence  Parameterization 


21.  Vertical  turbulence  is  handled  by  using  the  concept  of  eddy  viscos¬ 
ity  and  diffusivity  to  represent  the  velocity  and  density  correlation  terms 
that  arise  from  a  time  averaging  of  the  governing  equations.  These  eddy  coef¬ 
ficients  are  computed  from  mean  flow  characteristics  using  a  simplified 
second-order  closure  model  originally  developed  by  Donaldson  (1973).  The 
closure  model  has  been  further  developed  and  applied  to  various  types  of  flows 
by  Lewellen  (1977)  and  Sheng  (1982,  1986).  A  discussion  of  the  implementation 
of  the  turbulence  model  taken  from  Sheng*  follows. 

22.  Basically  all  second— order  closure  models  solve  the  Reynolds  stress 
equations  in  addition  to  the  equations  for  the  mean  flow.  These  equations  can 
be  found  in  any  textbook  on  turbulence  and  thus  are  not  reproduced  here. 
Assuming  local  equilibrium  of  turbulence,  the  Reynolds  stress  equations  can  be 
greatly  simplified  since  there  is  no  time  evolution  nor  spatial  diffusion  of 
the  second-order  correlations.  If  one  considers  only  the  vertical  turbulent 
transport,  these  equations  can  be  expanded  into 


*  Y.  P.  Sheng.  1990.  "A  Simplified  Second  Order  Closure  Model  of  Turbulent 
Transport."  Unpublished  paper  prepared  for  US  Army  Engineer  Waterways 
Experiment  Station,  Vicksburg,  MS. 
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^  {Fir  +  (2b-l)  ^  2  =  0  (15) 

^  ^  (2b-l)  2  II  =  0  (16) 

^  ^  (2b-l)  ^  4l  =  0  (17) 


where 


and 


(18) 
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where  u'  ,  v'  ,  and  w'  are  the  turbulent  velocity  fluctuations  and  the 
overbar  indicates  time  averaging.  These  equations  can  be  manipulated  to 
produce  an  equation  relating  the  turbulent  kinetic  energy  and  the  macroscale 
of  turbulence  A  to  the  mean  flow  shear  and  stratification  (given  by  the 
Richardson  number  Ri) ; 


SA^b^sQ'*  +  A[(bs  +  3b  +  7b2s)Ri  -  Abs(l  -  2b)]Q2 


(26) 


+  b(s  +  3  +  4bs)Ri^  +  (bs  -  A)(l  -  2b)Ri  =  0 
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23.  As  discussed  by  Sheng,*  it  can  also  be  shown  that  the  following 
relations  hold; 
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*  1990,  op.  cit. 
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where 


(j  = 


Ri 


and 
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Thus,  after  the  velocity  shear  and  flow  stratification  are  determined,  q  can 
be  computed  using  Equations  26  and  28.  w'  w'  is  then  determined  from 


_  ^ 

w'  w'  =  -n - 3- 

^(1  -  2w) 


(34) 


Finally,  after  A  is  prescribed,  u'  w'”  and  v'  w*  can  be  computed  from 
Equations  29  and  30  and  the  vertical  eddy  coefficients  can  be  determined  from 
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If  the  Richardson  Number  Ri  never  exceeds  a  critical  value  of  0.55  at  any 
vertical  position  in  the  water  column,  the  initial  distribution  of  A  at  a 
vertical  position  z  is  computed  from  0.65z  such  that  A  is  zero  at  both 
the  surface  and  the  bottom  with  its  maximum  value  occurring  at  middepth. 
However,  if  a  pycnocline  does  exist  (Ri  >  0.55),  then  the  initial  distribution 


of  A  is  computed  such  that  A  is  zero  at  the  surface,  at  the  pycnocline, 
and  at  the  bottom  with  maximum  values  occurring  halfway  between  the  surface 
and  the  pycnocline  and  the  bottom.  This  initial  distribution  is  then  modified 
by  the  three  basic  constraints  below: 


dA 
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(39) 


where  N  is  the  Brunt-Vaisala  frequency.  Equation  39  states  that  A  is  less 
than  a  fraction  of  the  spread  of  turbulence  as  measured  by  the  distance 
between  the  location  of  a  maximum  q^  to  where  q^  is  equal  to  25  percent  of 
the  maximum.  The  coefficient  Q<.ut  is  on  the  order  of  0.15  to  0.25. 
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PART  III:  APPLICATION  WITH  AN  EARLY  VERSION  OF  THE  MODEL  TO  THE 
JUNE-JULY  1980  AND  APRIL  1983  DATA  SETS 

Strategy 

24.  Verification  of  the  3— D  model  was  conducted  in  various  phases. 
Initially,  simplified  problems  for  which  closed  form  solutions  could  be  writ¬ 
ten  were  simulated  to  aid  in  assessing  model  performance  as  well  as  to  guide 
needed  model  modifications.  The  next  phase  involved  employing  data  that  had 
been  collected  on  a  physical  model  of  the  Chesapeake  Bay  in  Maryland.  How¬ 
ever,  because  these  data  do  not  contain  wind  effects,  the  effort  devoted  to 
this  phase  of  the  verification  process  was  minimal.  Ultimately  the  verifica¬ 
tion  of  any  numerical  model  must  be  conducted  using  sets  of  field  data.  In 
addition,  it  is  desirable  that  the  various  sets  of  field  data  represent  dif¬ 
ferent  physical  conditions  that  are  considered  important  and  routinely  occur 
in  the  system  being  modeled. 

25.  A  preliminary  study  of  synoptic  field  data  available  on  Chesapeake 
Bay  revealed  three  fairly  extensive  data  sets  that  were  each  approximately 

1  month  long.  As  described  in  paragraph  7,  the  first,  collected  during  June- 
July  1980  by  Dr.  Boicourt,  represents  dry  summer  conditions  on  Chesapeake  Bay. 
The  next  two  were  collected  by  the  NOS  during  1983.  One  set,  collected  during 
April  1983,  represents  a  relatively  large  freshwater  runoff  into  the  bay, 
whereas  the  other,  collected  during  September  1983,  captured  an  important 
wind-mixing  event.  These  data  sets  are  summarized  in  Table  1. 

26.  Model  development  during  this  study  was  an  ongoing  exercise,  e.g., 
the  rederivation  of  the  convective  terms,  the  incorporation  of  the  QUICKEST 
scheme,  the  incorporation  of  the  ability  to  compute  temperature,  and  the  solu¬ 
tion  of  the  governing  equations  on  the  z-plane  rather  than  the  sigma-plane. 
Therefore,  it  is  difficult  to  state  the  differences  between  the  version  (or 
versions)  of  the  model  applied  to  the  June-July  1980  and  April  1983  data  sets 
and  the  final  version  applied  to  the  September  1983  data  set  and  the  yearlong 
simulations  since  changes  were  made,  often  on  a  daily  basis.  However,  major 
differences  included  correcting  an  error  in  the  coding  of  the  advection  scheme 
near  inflow  boundaries,  applying  a  three-point  smoothing  equation  to  effec¬ 
tively  diffuse  the  vertical  turbulence,  computing  the  surface  heat  exchange 
coefficients  with  wind  data  adjusted  to  reflect  over— water  winds,  specifying 


20 


time-varying  salinity  and  temperature  at  the  ocean  boundary,  and  distributing 
ungaged  inflows  below  the  fall  lines  of  the  tributaries.  Since  these  individ¬ 
ual  differences  were  incorporated  over  a  period  of  time  during  many  simula¬ 
tions  of  the  September  1983  data  set,  it  is  difficult  to  ascertain  the  impact 
of  each  individually.  However,  an  inspection  of  the  September  1983  verifica¬ 
tion  results  reveals  that  collectively  they  resulted  in  a  much  better  matching 
of  the  computed  results  with  the  field  data. 

27.  Application  to  the  three  month-long  data  sets  are  discussed  in  two 
parts.  Results  from  the  June-July  1980  and  April  1983  applications  are  pre¬ 
sented  first  with  minimal  discussion  before  presenting  results  from  the 
September  1983  application.  Because  an  earlier  version  of  the  model  was 
applied  to  these  two  data  sets,  these  applications  should  not  be  viewed  as 
verification  efforts.  As  noted,  the  version  of  the  model  and  values  of  model 
coefficients  applied  to  the  September  1983  data  set  were  the  same  as  applied 
to  the  yearlong  simulations.  Thus,  this  application  is  considered  much  more 
important  and  as  a  result  is  discussed  in  more  detail  in  Part  IV.  Results 
from  the  June-July  1980  and  April  1983  data  sets  should  be  viewed  in  a 
qualitative  sense  as  demonstrating  basic  behavior  of  the  model  rather  than  a 
quantitative  verification  of  the  model. 

June-Julv  1980 

Description  of  field  data 

28.  As  discussed  by  Blumberg  et  al.  (1991),  velocity,  salinity,  and 
temperature  data  were  available  from  Boicourt  at  several  locations.*  Data  at 
the  four  stations  (M3,  WT4 ,  SP2 ,  BBl)  shown  in  Figure  3  were  selected  for 
comparison  with  model  results.  These  data  were  measured  at  either  a  10-min  or 
30-min  interval. 

29.  Wind  data  at  Norfolk  International  Airport  and  at  the  Baltimore- 
Washington  International  Airport  (BWl)  at  3— hour  intervals  were  available  for 
the  June-July  1980  period.**  The  wind  data  are  given  in  terms  of  wind 
direction  and  speed.  Hourly  tidal  height,  corrected  to  the  National  Geodetic 
Vertical  Datum  (NGVD) ,  were  available  at  six  locations:  Havre  de  Grace,  MD 


*  Personal  Communication  with  W.  Boicourt  of  Horn  Point  Environment  Labora¬ 
tories,  Cambridge,  MD. 

**  Personal  Communication  with  National  Climatic  Data  Center,  Asheville,  NC . 
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(HdG) ;  Colonial  Beach,  VA  (COL);  Annapolis,  MD  (ANNI ;  Solomons,  MD  (SOL); 
Hampton  Roads,  VA  (HRd) ;  and  Chesapeake  Bay  Tunnel,  VA  (CBT).*  The  locations 
of  these  gages  are  also  shown  in  Figure  3.  Average  daily  riverflow  data  at 
the  tributary  fall  lines  were  obtained  from  the  US  Geological  Survey  Water 
Resources  Data  reports  (USGS  1981a  and  b) . 

Initial  conditions 

30.  To  save  on  computation  time,  the  initial  3-D  salinity  and  tempera¬ 
ture  fields  were  constructed  to  be  as  close  to  reality  as  possible.  The  pro¬ 
cedure  used  is  discussed  in  more  detail  in  Part  IV.  The  initial  velocity 
field  was  taken  to  bf  zero  and  the  water  surface  was  taken  to  be  at  mean  sea 
level . 

Boundary  forcing  data 

31.  Boundary  forcings  consist  of  tides,  winds,  Freshwater  inflows,  and 
the  exchange  of  heat  at  the  water  surface.  The  tide  imposed  at  the  ocean 
boundary  (M3)  is  shown  in  Figure  4.  In  addition  to  the  tide  at  the  ocean 
boundary,  vertical  distributions  of  salinity  and  temperature  were  also  pre- 
toribed.  However,  in  this  application  these  distributions  were  not  allowed  to 
vary  with  time.  Thus,  on  inflow  the  values  prescribed  in  the  initial  condi¬ 
tions  became  the  boundary  conditions.  In  all  applications,  since  data  were 
not  available,  the  Chesapeake  and  Delaware  Canal  was  closed. 

2'*.  Daily  freshwater  ’nflows  (Figure  5)  were  prescribed  at  the  head  of 
all  major  tributaries.  All  of  these  flows  were  quite  low.  The  time-varying 
temperature  of  each  inflow  was  taken  to  be  the  equilibrium  temperature  com¬ 
puted  from  meteorological  data  at  the  Patuxent  River  Naval  Station  (PAX).** 

33.  Wind  stress  in  tne  lower  to  middle  bay  was  computed  from  linearly 
interpolated  wind  data  from  the  Norfolk  and  BWI  stations.  North  of  the  BWI 
station,  only  wind  data  from  the  BWI  station  were  used.  Wind  vectors  from  the 
Norfolk  and  BWI  stations  are  presented  in  Figures  6  and  7,  respectively. 

These  data  have  been  adjusted  to  reflect  over-water  winds.  This  was  accom¬ 
plished  by  multiplying  the  north-south  and  east-west  components  of  the  wind  at 
the  three  stations  by  factors  obt'>ined  from  Goodrich  (1985).  For  the  Norfolk 
winds  these  factors  were  1.53  and  1.85,  respectively;  and  for  the  Patuxent  and 
BWI  winds,  the  factors  used  were  2.03  and  1.40. 


*  Personal  Communication  with  National  Ocean  Service,  Rockville,  MD. 

**  Personal  Communication  with  National  Climatic  Data  Center,  Asheville,  NC . 
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34.  As  previously  noted,  the  surface  heat  exchange  is  handled  through 
the  concept  of  an  equilibrium  temperature  and  a  surface  heat  exchange  coeffi¬ 
cient  as  discussed  by  Edinger,  Brady,  and  Geyer  (1974).  These  were  computed 
from  meteorological  data  taken  at  the  Patuxent  Naval  Air  Station  and  are 
presentevj  in  Table  2.* 

Verification  results 

35.  With  the  boundary  forcings  and  initial  state  described  previously, 
computations  were  made  on  the  numerical  grid  shown  in  Figure  2.  The  initial 
salinity  and  temperature  fields  were  frozen  for  the  first  7  days.  Computed 
salinity,  temperature,  and  velocity  were  saved  at'  the  locations  previously 
noted.  Computed  water-surface  elevations  were  saved  at  the  interior  tide 
stations  noted  in  Figure  3. 

36.  Comparisons  of  the  computed  and  recorded  tide  at  Hampton  Roads,  VA, 
and  Solomons,  Annapolis,  and  Havre  de  Grace,  MD,  are  presented  in  Figures  8- 
11.  Generally,  the  amplitudes  match  well  and  reflect  wind  effects  as  well  as 
the  ast’'r>'' ^.iiical  tide.  Phasing  differences  can  be  seen  in  the  results  at 
Havre  de  Grace  (Figure  11). 

37.  A  comparison  at  three  depths  of  the  major  velocity  component  at  the 
ocean  entrance  (M3)  is  presented  in  Figure  12.  A  similar  plot  at  three  depths 
of  the  along-bay  component  of  the  flow  velocity  at  station  WT4  is  given  in 
Figure  13.  In  these  plots,  the  u-velocity  is  he  along-bay  component,  and  the 
v-velocity  is  the  across-bay  component.  Comparisons  of  computed  and  recorded 
salinities  and  temperature  at  station  WT4  are  given  in  Figures  14  and  15, 
respectively.  Comparisons  of  computed  and  recorded  velocity,  salinity,  and 
temperature  at  station  SP2  are  presented  in  Figures  16-18,  and  at  station  BBl 
in  Figures  19-21 . 

38.  Velocity  comparisons  are  generally  favorable,  especially  if  one 
considers  that  the  recorded  velocity  occurs  at  a  point,  whereas  the  model 
velocity  represents  an  average  over  a  large  computational  cell.  The  salinity 
and  tempera'^ure  comparisons  are  not  as  good.  However,  as  previously  noted, 
these  results  were  generated  with  an  early  version  of  the  model.  The  excel¬ 
lent  results  presented  in  Part  IV  for  the  September  1983  data  set  are  from  the 
latest  version  of  the  model  in  which  the  various  modifications  discussed  in 
paragraph  26  have  been  made.  Once  again  it  should  be  noted  that  the  surface 
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heat  exchange  data  presented  in  Table  1  were  gene''ated  without  correcting  the 
Patuxent  wind  data  to  account  for  winds  over  water.  This  correction  resulted 
in  a  much  better  comparison  of  temperatures  during  the  September  1983  applica¬ 
tion.  Due  to  time  and  funding  constraints,  the  1980  data  set  was  not  rerun 
using  the  latest  version  of  the  model. 


April  1983 

39.  Data  for  this  application  were  obtained  from  a  circulation  survey 
of  Chesapeake  Bay  conducted  by  NOS  during  April  1983*  and  are  discussed  in 
detail  by  Blumberg  et  al.  (1991).  This  data  set  represents  conditions  during 
a  large  spring  runoff  ev'ent . 

Description  of  field  data 

40.  As  shown  in  Figure  22,  current  meter  data  were  available  for  four 
stations;  bay  entrance  (BE),  midbay  (MB),  Wolf  Trap  (WT) ,  and  Bay  Bridge 
(BB) .  These  data  include  current  direction  and  amplitude,  water  temperature, 
pressure,  and  conductivity.  The  data  were  measured  at  a  10-min  interval.  The 
temperature  and  conductivity  were  used  to  compute  salinity  from  the  following 
expression**; 

S  -=  (([(-1.02527  *  RF  +  4.81236)  *  RF  -  9.04061] 

*  RF  +  11.95364)  *  RF  +  28.29988)  *  RF  (40) 

where 

RF  =  RT  +  RT  *  RD  *  TF  *  (BR  +  TF  +  FRT)  *  10'^ 

and 

RT  =  C/CKT 
RD  =  RT  -  1 .0 

BR  =  ([(-26.9  *  RD  +  3.09)  *  RD  -  8.52]  *  RD  +  67.1) 

FRT  =  -0.25  *  RD  *  TF 

with 

CKT  -  ([(-2.217  *  10  *  T  -  2.5813  *  10“^) 

*  T  +  4.6704  *  10*3]  *  X  +  0.86062)  *  T  +  29.0473 


*  Personal  Communication  with  National  Ocean  Service,  Asheville,  NC . 
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TF  =•  T  -  15.0 


and 

T  =  temperature,  °C 
C  “  conductivity,  mmhos/cm 
S  =  salinity,  ppt 

41.  Hourly  wind  data  at  the  Patuxent  River  Naval  Air  Station  and  the 
Norfolk  International  Airport  were  available  for  the  April  1983  period.* 

These  data  include  wind  direction  and  wind  amplitude.  Figures  23  and  24  pres¬ 
ent  wind  data  from  these  two  stations.  Locations  of  the  stations  are  shown  in 
Figure  22.  These  data  have  been  adjusted  to  represent  over-water  winds  as 
discussed  in  paragraph  33. 

42.  Daily  riverflow  data  at  the  tributary  fall  lines  were  obtained  from 
the  US  Geological  Survey  Water  Resources  Data  reports  for  the  April  1983 
period  (USGS  1984a  and  b) .  These  data  are  presented  in  Figure  25  for  the 
James,  York,  Rappahannock,  Potomac,  Patuxent,  Patapsco,  Susquehanna,  and 
Choptank  Rivers. 

43.  Hourly  tidal  height  data  at  six  locations  were  available  from  NOS 
for  the  April  1983  period.**  The  location  of  these  gages  is  given  in  Fig¬ 
ure  22.  The  values  have  been  corrected  to  NGVD.  The  boundary  tide  at  CBT  is 
given  in  Figure  26. 

44.  Equilibrium  temperatures  and  surface  heat  transfer  coefficients 
computed  from  meteorological  data  at  the  Patuxent  River  Naval  Air  Station  for 
the  April  1983  period  are  presented  in  Table  3.*  As  with  the  1980  data  set, 
the  Patuxent  winds  were  not  adjusted  to  reflect  over-water  winds  in  the 
surface  heat  exchange  computations. 

Verification  results 

45.  Initial  conditions  were  generated  as  in  the  1980  application 
(paragraph  30) . 

46.  Comparisons  of  the  recorded  and  computed  tides  at  Hampton  Roads  and 
Colonial  Beach,  VA,  and  Solomons,  Annapolis,  and  Havre  de  Grace,  MD ,  are  pre¬ 
sented  in  Figures  27-31,  respectively.  Results  are  quite  good  except  at 


*  Personal  Communication  with  National  Climatic  Data  Center,  Rockville,  MD. 

**  Personal  Communication  with  National  Ocean  Service,  Asheville,  NC . 
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Havre  de  Grace.  The  reason  for  the  much  larger  computed  tidal  range  at 
Havre  de  Grace  over  the  last  half  of  the  month  is  not  clear.  This  behavior  at 
Havre  de  Grace  was  not  nearly  as  pronounced  in  the  1980  application 
(Figure  11) . 

47.  Comparisons  of  computed  and  recorded  velocities  at  stations  BE  and 
WT  are  presented  in  Figures  32  and  33,  and  comparisons  of  velocity,  salinity, 
and  temperature  at  stations  MB  and  BB  (Figure  22)  are  presented  in  Figures  34- 
39.  Again,  it  should  be  remembered  that  recorded  velocities  are  point  values, 
whereas  computed  velocities  are  all  cell— averaged  values.  In  addition,  the 
previous  comments  concerning  the  use  of  an  early  version  of  the  model  and 
erroneous  surface  heat  exchange  data  are  applicable. 

Conclusions 


48.  Based  upon  an  overall  comparison  of  model  results  with  field  data 
at  several  locations,  the  general  conclusion  after  these  two  applications  was 
that  the  3-D  hydrodynamic  model  responded  to  boundary  forcings  in  a  reasonable 
fashion.  After  correcting  model  leficiencies  during  the  September  1983  appli¬ 
cation  (Part  IV)  and  after  accounting  for  ungaged  tributary  inflows  and  the 
adjustment  of  the  Patuxent  wind  data  in  the  computation  of  the  surface  heat 
exchange,  comparisons  of  computed  and  field  data  improved  substantially  for 
the  September  1983  application.  As  a  result,  confidence  in  the  ability  of  the 
model  to  accurately  compute  the  3-D  hydrodynamics  of  the  Chesapeake  Bay  was 
increased.  This  confidence  was  increased  even  more  during  the  yearlong 
simulations  discussed  in  Part  V. 
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PART  IV:  VERIFICATION  OF  THE  FINAL  MODEL  TO  THE 
SEPTEMBER  1983  DATA  SET 

49.  These  data  illustrate  the  destratification  of  the  Chesapeake  Bay 
waters  during  a  strong  wind  event  with  a  subsequent  restratification  as  the 
wind  abated.  In  addition,  rapid  cooling  of  the  surface  waters  near  the  middle 
of  the  month  contributed  to  the  destratification  process. 

Description  of  Field  Data 

50.  As  with  the  April  1983  data,  current  meter  data  were  available  for 
four  long-term  stations:  bay  entrance  (BE),  Woif  Trap  (WT) ,  midbay  (MB),  and 
Bay  Bridge  (BB).*  These  locations  arc  shown  in  Figure  40.  These  data  in¬ 
clude  current  direction  and  amplitude,  water  temperature,  pressure,  and  con¬ 
ductivity.  They  were  measured  at  a  10-min  time  interval.  The  temperature  and 
conductivity  were  used  to  compute  salinity  from  Equation  40. 

51.  Hourly  wind  data  at  the  Patuxent  River  Naval  Air  Station  and  wind 
data  at  Norfolk  International  Airport  at  3-hour  intervals  were  available  for 
the  September  1983  period.**  The  wind  data  are  given  in  terms  of  wind 
direction  and  speed.  Hourly  tidal  height,  corrected  to  NGVD,  at  seven  loca¬ 
tions  were  obtained  from  NOS:  Havre  de  Grace,  MD  (HdG) ,  Colonial  Beach,  VA 
(COL),  Lewisetta,  VA  (LEW),  Annapolis,  MD  (ANN) ,  Solomons,  MD  (SOL),  Hampton 
Roads,  VA  (HRd) ,  and  Chesapeake  Bay  Tunnel,  VA  (CBT).*  The  locations  of  these 
gages  are  also  shown  in  Figure  40.  Average  daily  riverflow  data  at  the  tribu¬ 
tary  fall  lines  were  obtained  from  the  US  Geological  Survey  Water  Resources 
Data  reports  (USGS  1984a  and  b) . 

Initial  Conditions 

52.  To  save  on  computation  time,  the  initial  3-D  salinity  and  tempera¬ 
ture  fields  were  constructed  to  be  as  close  to  reality  as  possible.  These 
were  established  by  using  the  available  field  data.  The  value  for  each  indi¬ 
vidual  cell  of  the  3-D  grid  was  first  set  to  be  that  of  the  nearest  field  data 
point,.  The  resulting  3-D  fields  of  salinity  and  temperature  were  then 
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smoothed  several  times  using  a  three— point  smoothing  equation.  Figures  41  and 
42  show  the  vertical  structure  of  the  prescribed  salinity  field  at  stations  MB 
and  BB.  Based  upon  the  limited  field  data  available,  it  can  be  seen  that  ini¬ 
tial  stratification  of  the  bay  is  represented  well.  The  initial  velocity 
field  was  taken  to  be  zero  and  the  water  surface  was  taken  to  be  at  mean  sea 
level . 


Boundary  Forcing  Data 

53.  As  for  the  1980  and  April  1983  applications,  boundary  forcings 
consist  of  tides,  winds,  freshwater  inflows,  and  the  exchange  of  heat  at  the 
water  surface.  The  tide  imposed  at  the  ocean  boundary  (CBT)  is  shown  in  Fig¬ 
ure  43.  In  addition  to  the  tide  at  the  ocean  boundary,  time— varying  vertical 
distributions  of  salinity  and  temperature  were  also  prescribed  using  the  field 
data  at  station  BE.  Recall  that  in  the  previous  applications  these  vertical 
distributions  were  constant  in  time.  As  previously  noted,  since  data  were  not 
available,  the  Chesapeake  and  Delaware  Canal  was  closed. 

54.  Daily  freshwater  inflows  were  prescribed  at  the  head  of  all  major 
tributaries.  All  of  these  flows  were  quite  low,  with  over  half  of  the  total 
freshwater  inflow  coming  from  the  Susquehanna  River.  These  inflows  are  shown 
in  Figure  44.  The  time-varying  temperatures  of  these  inflows  were  also 
prescribed.  The  salinity  of  the  tributary  inflows  was  taken  to  be  zero. 

55.  Wind  stress  in  the  lower  to  middle  bay  was  computed  from  linearly 
interpolated  wind  data  from  the  Norfolk  and  Patuxent  stations.  North  of  the 
Patuxent  station,  only  wind  data  from  the  Patuxent  station  were  used.  Wind 
vectors  from  these  two  stations  are  presented  in  Figures  45  and  46  and  have 
been  adjusted  to  represent  over-water  winds. 

56.  As  previously  noted,  the  surface  heat  exchange  was  computed  as 
described  in  paragraph  34.  Unlike  the  previous  two  applications  discussed, 
the  wind  data  employed  in  the  computation  of  the  surface  heat  exchange  data 
presented  in  Table  4  were  adjusted  to  represent  over— water  winds. 

Model  Coefficients 


57.  The  horizontal  eddy  viscosity  and  diffusivity  coefficients  are 
prescribed  along  with  parameters  connected  with  the  computation  of  the 
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vertical  turbulence  coefficients  and  the  bottom  drag  coefficient.  The  bottom 
roughness  was  set  such  that  the  bottom  drag  coefficient  was  computed  to  have  a 
value  of  0.0028  throughout  the  bay.  Values  for  the  horizontal  eddy  viscosity 
and  diffusivity  were  taken  to  be  zero.  Background  values  for  the  vertical 
eddy  viscosity  and  diffusivity  were  set  to  be  10  cmVsec  and  0.005  cm^/sec, 
respectively.  These  same  values  were  used  in  the  previous  two  applications. 

5b.  As  discussed  in  Paj.t  II,  the  vertical  Luibulence  model  contains  one 
free  parameter,  that  can  be  varied.  However,  its  value  should  be  in  the 

range  of  0.15—0.25.  In  this  range,  the  computed  eddy  coefficients  were  fairly 
insensitive  to  variations  in  •  Thus,  was  set  to  0.20.  The  other 

three  parameters,  b  ,  A  ,  and  S  ,  are  considered  to  be  constants.  However, 
several  runs  with  the  September  1983  data  set  were  made  in  which  the  values  of 
these  parameters  were  varied.  The  basic  conclusion  drawn  from  these  results 
was  that  the  default  values  should  not  be  changed  and  these  parameters  should 
indeed  be  considered  model  constants. 


Verification  Results 


59.  With  these  boundary  forcings  and  initial  state,  computations  were 
made  on  the  numerical  grid  shown  in  Figure  2.  The  initial  salinity  and  tem¬ 
perature  fields  were  frozen  for  the  first  5  days.  Computed  salinity,  tempera¬ 
ture,  and  velocity  were  saved  at  the  locations  of  the  current  meter  stations 
in  Figure  40  for  comparison  with  the  field  data.  Computed  water— surface  ele¬ 
vations  were  saved  at  the  six  interior  tide  stations  noted  in  Figure  40. 
Furthermore,  monthly  averaged  velocity  vectors  were  computed  for  comparison 
with  observed  long-term  circulation  patterns  in  Chesapeake  Bay. 

60.  In  addition  to  time  series  plots  showing  a  visual  comparison  of 
model  results  to  field  data,  difference  measures  have  also  been  computed. 

Bias  is  described  by  the  mean  bias  error  (MBE)  defined  as 


where 


MBE  - 


E  (Mi  -  Oi) 

1.1 

N 


(91) 


N  ”  number  of  data  points 
Mi  -  model  results 
0.  -  observed  data 
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Average  difference  is  described  by  the  root  mean  square  error  (RMSE)  or 
perhaps  the  mean  absolute  error  (MAE).  These  are  defined  by  the  following: 


RMSE 


E  (Mi  -  Oi)" 

- R - 


(42) 


MAE 


H 


E  K  -  Oi 


N 


(43) 


Finally,  the  average  relative  error  (ARE)  is  defined  as 


ARE 


E  |Mi  -  Oil 

- 

E  Oi 

i-l 


(44) 


61.  Comparison  of  the  computed  and  recorded  tides  at  Hampton  Roads, 
Lewisetta,  and  Colonial  Beach,  VA,  and  Solomons,  Annapolis,  and  Havre  de 
Grace,  MD,  are  presented  in  Figures  47—52.  Table  5  presents  a  comparison  of 
model  and  observed  maximum  flood  and  minimum  ebb  water  surface  elevations, 
along  with  a  comparison  of  the  timing  of  the  peaks,  at  these  stations.  Low 
water  at  Hampton  Roads  is  consistently  computed  to  be  too  high;  however,  the 
range  at  most  other  stations  is  about  right,  except  for  Havre  de  Grace.  As 
with  all  the  simulations,  the  computed  range  at  that  station  is  too  high. 
Around  the  middle  of  the  month  the  computed  water  surface  at  most  stations 
experiences  too  much  setup.  This  is  probably  due  to  an  inaccuracy  in  the  wind 
data  or  perhaps  a  barometric  disturbance  that  is  not  modeled.  The  effect  of 
the  wind  shift  around  the  20th  of  the  month  is  quite  evident  in  the  water 
surface  plots  in  the  upper  bay.  Figures  53—58  are  similar  water— surface  plots 
but  with  all  frequencies  below  36  hours  filtered  from  the  results.  Thus, 
these  plots  illustrate  a  comparison  of  recorded  and  computed  water-surface 
elevations  due  to  nonastronomical  events,  e.g.,  wind  and  freshwater  runoff. 

The  difference  measures  for  the  surface  elevations  are  presented  in  Table  6. 
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The  bias  is  positive  at  all  stations  and  is  generally  about  0.05  m  or  less. 

The  difference  error  is  less  than  about  0.1  m  at  all  stations  except  Havre  de 
Grace.  This  is  a  unique  station  located  at  the  head  of  the  bay  in  the 
Susquehanna  Flats  which  are  not  represented  well  in  the  model. 

62.  Figures  59—62  illustrate  the  ability  of  the  numerical  model  to 
reproduce  flow  velocity  well.  The  impact  of  the  sudden  shift  in  the  wind 
around  the  20th  of  the  month  can  clearly  be  seen  at  stations  MB  and  BB. 

Table  7  presents  a  comparison  of  model  and  observed  along— bay  currents  near 
the  surface  at  maximum  flood  and  minimum  ebb.  Table  8  presents  the  difference 
measures  for  the  velocity.  From  an  inspection  of  the  time  series  plots,  it 
appears  the  difference  error  is  caused  more  by  phase  errors  than  by  range 
errors . 

63.  Salinity  and  temperature  results  are  presented  at  the  three 
interior  stations  shown  in  Figure  40  in  Figures  63—65  and  66-68,  respectively. 
The  salinity  plots  at  stations  MB  and  BB  demonstrate  that  the  vertical  turbu¬ 
lence  closure  model  responds  reasonably  well  to  the  wind-mixing  event.  As  a 
result  of  the  turbulence  generated  by  the  velocity  shear  created  by  the  wind, 
the  bay  destratified  but  then  began  to  restratified  within  a  few  days.  To 
adequately  capture  this  event,  it  was  found  that  some  diffusion  of  the  turbu¬ 
lence  generated  was  required.  This  was  accomplished  by  using  a  three— point 
smoothing  equation  on  the  eddy  viscosity  and  diffusivity  coefficients 
generated  by  the  turbulence  model  discussed  in  Part  II.  With  k  being  the 
vertical  layer  indicator  and  D  representing  either  the  eddy  viscosity  or 
diffusivity,  the  following  equation  is  employed: 


D.  = 


Dk-i  ^ 

3 - 


(45) 


Figure  69  shows  the  salinity  at  the  MB  station  without  the  three-point  smooth¬ 
ing.  A  comparison  with  Figure  64b  clearly  demonstrates  the  impact  of  employ¬ 
ing  such  a  smoothing  scheme.  Difference  measures  for  the  salinity  are 
presented  in  Table  9.  Generally,  the  mean  bias  is  less  than  1  ppt  with  the 
difference  error  also  less  than  1  ppt,  except  at  Wolf  Trap.  Difference 
measures  for  the  temperature  are  presented  in  Table  10.  Generally,  the 
absolute  error  as  well  as  the  mean  bias  error  is  less  than  0.5°C. 


31 


64.  Since  flow  fields  generated  by  the  hydrodynamic  model  are  to  be 
used  in  the  long-term  computation  of  water  quality  parameters,  it  is  important 
to  demonstrate  the  ability  of  the  model  to  compute  the  proper  residual  circu¬ 
lation  of  the  bay.  Figure  70  shows  the  computed  monthly  averaged  near-surface 
and  near- bottom  currents  in  the  bay.  Note  the  classical  gravitational  circu¬ 
lation,  with  the  surface  water  flowing  toward  the  ocean  and  bottom  water  mov¬ 
ing  up  the  bay.  Figure  71  is  a  similar  plot  that  was  constructed  by  Blumberg, 
from  records  of  field  data  collected  by  NOS  at  various  locations  over  a  period 
of  time  from  1976  to  1983.*  Only  those  records  of  at  least  15  days  in  length 
were  employed  in  constructing  these  plots.  Note  that  the  magnitude  of  the 
model-computed  near-surface  and  near-bottom  residual  currents  is  in  the  5-  to 
10-cm/sec  range  as  is  that  for  the  field  data.  Although  one  should  not  neces¬ 
sarily  expect  results  representing  an  average  of  data  collected  from  1976  to 
1983  to  agree  with  the  computed  monthly  averaged  circulation  for  September 
1983,  it  can  be  seen  that  in  a  qualitative  sense  these  results  do  justify 
placing  confidence  in  the  model's  ability  to  compute  the  residual  circulation 
of  the  bay. 

65.  Computed  nondimensional  eddy  diffusivity  at  the  midbay  and  Bay 
Bridge  stations  before  the  use  of  Equation  45  are  presented  in  Figures  72  and 
73.  The  values  plotted  should  be  multiplied  by  10  cmVsec  to  yield  dimen¬ 
sional  values.  An  inspection  of  the  computed  coefficients  shows  values  of 
perhaps  500-1,000  cmVsec  being  sustained  for  a  couple  of  days  during  the 
wind-mixing  event.  These  plots  also  illustrate  the  importance  of  shear 
stresses  at  the  boundaries  in  the  generation  of  turbulence. 

Conclusions 


66.  Based  upon  an  overall  comparison  of  model  results  with  available 
field  data  at  several  locations,  the  general  conclusion  is  that  the  3-D  model 
is  a  good  representation  of  the  hydrodynamics  of  the  Chesapeake  Bay.  This  is 
further  qualitatively  substantiated  by  the  excellent  agreement  of  the  computed 
monthly  averaged  flow  field  with  average  circulation  values  determined  from 
field  data  collected  by  NOS  during  1976-83. 


*  Personal  Communication  with  Dr.  Alan  Blumberg,  HydroQual ,  Inc.,  Mahwah ,  NJ . 
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PART  V:  MODEL  VERIFICATION  ON  AN  ANNUAL  TIME  SCALE 


67.  Hydrodynamics  for  the  years  of  1984,  1985,  and  1986  are  required  as 
input  to  the  3-D  water  quality  model  of  Chesapeake  Bay.  These  years  represent 
high,  low,  and  average  freshwater  inflow  years,  respectively.  A  1— year  simu¬ 
lation  of  the  hydrodynamics  using  a  computational  time-step  of  5  min  takes 
about  10  hr  on  a  CRAY  Y— MP  computer.  Sequences  of  these  hydrodynamics  will  be 
constructed  to  drive  long-term  (decades)  water  quality  simulations. 

68.  For  each  of  these  years,  near-surface  and  near-bottom  salinity  and 
temperature  data  at  approximately  2-week  intervals  were  available  throughout 
the  bay  and  ics  major  tributaries.*  These  data,  as  well  as  tide  data,** 
were  used  as  additional  verification  of  the  hydrodynamic  model. 

69.  Time  series  plots  of  water  level,  salinity,  and  temperature  as  well 
as  seasonally  averaged  longitudinal  transects  of  salinity  along  the  main  bay 
and  its  major  tributaries  have  been  constructed.  The  location  of  these 
transects  is  given  in  Figure  74. 


1984  Simulation 


Boundary  conditions 

70.  The  1984  year  was  broken  into  five  seasons  as  follows  for  the  pur¬ 
pose  of  generating  seasonally  averaged  longitudinal  transects  of  salinity; 

Season  1  =»  1  Jan-25  Mar 
Season  2  =>  26  Mar-10  Jun 
Season  3  =>  11  Jun-27  Aug 
Season  4  =>  28  Aug— 25  Nov 
Season  5  =»  26  Nov— 31  Dec 

Time  series  results  and  boundary  forcing  data  are  presented  for  a  portion  of 
the  year.  The  ocean  boundary  tide  at  the  Chesapeake  Bay  Tunnel  is  shown  in 
Figure  75.  Wind  forcing  data  that  have  been  corrected  to  reflect  over-water 
winds  based  on  recorded  winds  at  the  Norfolk  and  BWI  Airports  are  presented  in 


*  Personal  Communication,  1989,  with  Chesapeake  Bay  Liaison  Office,  US  Envi¬ 
ronmental  Protection  Agency,  Annapolis,  MD. 

**  Personal  Communication  with  National  Ocean  Service,  Asheville,  NC. 
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Figures  76  and  77.*  Freshwater  inflows  on  the  James,  Potomac,  Patuxent,  and 
Susquehanna  Rivers  are  shown  in  Figure  78  (USGS  1985a  and  b) .  Surface  heat 
exchange  data  for  season  5  are  listed  in  Table  11. 

Results 

71.  Locations  for  which  results  are  presented  are  shown  in  Figure  79. 
Complete  results  are  given  in  Appendix  A.  A  typical  comparison  of  computed 
and  recorded  water  level  at  Annapolis,  MD,  is  given  in  Figure  80.  Figure  81 
is  an  expanded  view  of  this  comparison  near  day  301  for  two  tidal  cycles. 

72.  Salinity  comparisons  are  given  at  four  main  bay  stations  as  well  as 
stations  on  the  James,  Potomac,  and  Patuxent  Rivers  for  the  last  3  months  of 
1984.  Station  EE  3.5  was  located  in  a  relatively  shallow  area,  about  4.5  m, 
on  the  eastern  side  of  the  bay.  As  can  be  seen  in  Figure  82,  very  little 
stratification  ever  appears  in  either  the  field  data  or  model  results.  The 
Wolf  Trap  station  (station  CB  6.3)  was  in  the  same  location  as  in  the 
September  1983  data  set.  Figure  83  indicates  the  water  column  can  be  alter¬ 
nately  well  mixed  and  partially  stratified.  The  raidbay  station,  CB  5.1,  was 
located  in  the  deep  natural  channel  traversing  Chesapeake  Bay  and  thus  basi¬ 
cally  remained  stratified  throughout  the  year,  although  the  degree  of  strati¬ 
fication  depended  upon  the  degree  of  turbulence  and  the  amount  of  freshwater 
inflow.  Figures  84  and  85  show  the  ability  of  the  3-D  model  to  reproduce  the 
stratification  throughout  the  year.  The  bay  bridge  station,  CB  3.3W,  was 
located  nearer  the  saline  front  and  thus  experienced  greater  fluctuation.  As 
illustrated  in  Figure  86,  the  water  column  at  this  station  experiences  a  great 
deal  of  apparent  mixing  and  subsequent  restratification.  The  numerical  model 
reproduces  these  events  extremely  well.  Difference  measures  for  the  main  bay 
stations  are  presented  in  Table  12.  The  average  relative  error  is  about 

10  percent  for  the  salinity  computations. 

73.  Figure  87  illustrates  the  comparison  of  computed  and  recorded 
salinity  near  the  mouth  of  the  James  River.  Model  results  indicate  the  water 
column  is  stratified  throughout  the  year.  However,  some  of  the  field  data 
imply  a  well-mixed  condition  at  times.  The  exact  reason  for  this  discrepancy 
is  unknown.  The  coarseness  of  the  grid  prohibits  modeling  localized  topo¬ 
graphic  effects.  In  addition,  errors  in  the  specification  of  the  local  wind 
may  be  a  factor. 


*  Personal  Communication  with  National  Climatic  Data  Center,  Rockville,  MD . 
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74.  Figures  88  and  89  show  salinity  comparisons  on  the  Potomac  ( LE  2.2^ 
and  Patuxent  (LE  1.1)  Rivers,  respectively.  These  results  serve  to  demon¬ 
strate  that  conditions  on  the  tributaries  are  represented  reasonably  well .  An 
inspection  of  the  numerical  grid  in  Figure  2  reveals  the  rather  coarse  discre¬ 
tization  of  the  tributaries.  With  this  grid  resolution  it  is  encouraging  that 
conditions  are  reproduced  as  well  as  they  are  on  the  tributaries.  The  spa¬ 
tially  third-order  QUICKEST  scheme  employed  in  the  salt  and  temperature 
transport  equation  is  primarily  responsible  for  this  success.  Difference 
measures  on  the  tributaries  are  given  in  Table  12  and  show  that  the  average 
relative  error  is  approximately  15-20  percent. 

75.  Temperature  comparisons  at  the  same  stations  are  given  in  Fig¬ 
ures  90-96.  These  results  demonstrate  the  adequacy  of  the  surface  heat  ex¬ 
change  formulation  in  modeling  temperatures  in  Chesapeake  Bay.  With  the  bay 
being  so  large,  if  more  than  one  meteorological  station  had  been  employed, 
these  results  would  likely  have  been  even  better.  Difference  measures  for 
both  the  main  bay  and  the  tributaries  are  given  in  Table  13.  It  can  be  seen 
that  the  average  relative  error  is  generally  less  than  10  percent  in  the  main 
bay  and  less  than  15  percent  in  the  tributaries. 

76.  Seasonally  averaged  longitudinal  transects  of  salinity  have  been 
constructed  from  both  model  results  and  the  field  data.  The  location  of  these 
transects  are  shown  in  Figure  74.  Results  for  seasons  4  and  5  for  both  main 
bay  transects  and  the  James  and  Potomac  Rivers  are  given  in  Figures  97-100. 
When  viewing  these  results  it  should  be  remembered  that  the  seasonally  aver¬ 
aged  field  data  generally  involve  only  four  or  five  values  and  in  some  cases 
only  one  value.  The  range  of  the  field  data  is  shown  by  the  vertical  lines 
drawn  through  tlie  average  value  at  each  location  where  data  were  available. 

The  locations  of  the  mouths  f  the  Potomac  and  Patapsco  Rivers  are  noted. 

77.  The  final  set  of  plots  of  1984  results  are  presented  in  Fig¬ 
ures  101-104.  Data  from  all  of  the  main  bay  stations  and  all  of  the  tributary 
stations  have  been  grouped  together  separately  to  illustrate  a  comparison  of 
computed  and  recoided  stratification,  as  reflected  by  the  difference  in  salin¬ 
ity  between  the  top  and  bottom  layers,  for  the  complete  year.  As  shown  in 
Figures  101  and  103,  a  linear  ^egression  analysis  shows  the  square  of  the 
correlation  coefficient  R‘  to  be  above  0.80  for  both  the  m.ain  bay  and  tribu¬ 
tary  stations.  An  inspection  of  Figures  102  and  104  reveals  that  the  error  in 
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the  computed  stratification  is  less  than  about  1  ppt  for  approximately  70  per¬ 
cent  of  the  time. 


1985  Simulation 


Boundary  conditions 

78.  The  1985  year  was  broken  into  the  following  five  seasons  for  the 
purpose  of  generating  seasonally  averaged  longitudinal  transects  ot  salinity: 

Season  1  =»  1  Jan-28  Feb 
Season  2  =>  1  Mar-30  Apr 
Season  3  =»  1  May-15  Jul 
Season  4  =>  16  Jul— 18  Sep 
Season  5  =>  19  Sep-31  Dec 

79.  Boundary  forcing  data  are  presented  for  season  4  of  1985.  The 
ocean  boundary  tide  at  the  Chesapeake  Bay  Tunnel  is  given  in  Figure  105  with 
the  wind  data  at  the  Norfolk  and  BWI  airports  presented  in  Figures  106  and 
107.  Freshwater  inflows  for  season  4  on  the  James,  Potomac,  Patuxent  and 
Susquehanna  Rivers  are  shown  in  Figure  108  (USGS  1986a, b) .  Surface  heat  ex¬ 
change  data  for  season  4  are  listed  in  Table  14. 

Results 

80.  Locations  for  which  results  are  presented  are  shown  in  Figure  109. 
Complete  results  are  given  in  Appendix  B.  A  comparison  of  computed  and 
recorded  water  level  at  Colonial  Beach,  VA,  is  given  in  Figure  110,  This 
comparison  is  typical  of  results  at  other  stations.  Figure  111  is  an  expanded 
view  of  this  comparison  near  day  201  for  two  tidal  cycles. 

81.  Salinity  comparisons  are  given  at  four  main  bay  station  5  as  well  as 
stations  on  the  James,  Potomac,  and  Patuxent  Rivers  for  season  4.  These  sta¬ 
tions  are  the  same  ones  for  which  1984  results  were  presented.  Plots  showing 
salinity  comparisons  at  the  main  bay  stations  are  given  in  Figures  112-115.  A 
comparison  of  Figure  114  and  Figure  84  reveals  a  greater  stratification  in  the 
midbay  during  1984.  This  is  because  1984  was  a  higher  freshwater  inflow  year. 
Note  that  the  numerical  model  computes  the  reduced  stratification  in  1985 
quite  well.  The  apparent  mixing  and  subsequent  restratification  observed  in 
the  1984  data  at  the  Bay  Bridge  station  can  also  be  seen  in  the  1985  results 
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presented  in  Figure  115.  Cnee  again  the  model  does  an  excellent  job  of 
reproducing  these  events. 

82.  Salinity  competrisons  on  the  James,  Potomac,  and  Patuxent  Rivers  are 
presented  in  Figures  116-118.  Generally,  these  results  show  a  good  reproduc¬ 
tion  of  conditions  on  the  James  and  Patuxent  Rivers.  "t  does  appear  from 
Figure  117  that  higher  saline  conditions  than  are  reflected  in  the  field  data 
are  computed  in  the  lower  Potomac  River  during  season  4.  An  inspection  of 
results  at  tuis  station  for  the  complete  year  (see  Appendix  B)  reveals  that 
this  is  true  for  most  of  the  year.  The  coarseness  of  the  grid  in  the  tribu¬ 
taries  which  results  in  an  inadequate  representation  of  the  channel  is  proba¬ 
bly  responsible  for  the  increased  error  at  some  tributary  locations.  Tempera¬ 
ture  results  at  the  same  stations  are  presented  in  Figures  119-125  for  season 
4.  Difference  measures  for  both  salinity  and  temperature  are  given  in 
Tables  12  and  13.  It  appears  the  model  did  a  slightly  better  job  in  1984. 
However,  it  should  be  noted  that  fewer  field  data  were  available  in  1984  for 
comparison . 

83.  During  1985,  salinity  and  temperature  data  were  collected  at  1-m 
intervals  at  several  locations.*  These  have  been  used  along  with  model 
results  to  construct  Figure  126  showing  a  comparison  of  computed  and  recorded 
vertical  profiles  of  water  density  at  station  CB  5.1  at  six  times  during  the 
year.  These  results  demonstrate  that  the  model  maintains  the  proper  vertical 
stratification  with  the  pycnocline  in  approximately  the  correct  location. 

84.  Seasonally  averaged  salinities  have  been  computed  along  the 
tiansects  shown  in  Figure  74.  Results  for  season  4  for  both  main  bay 
transects  as  well  as  the  James  and  Potomac  Rivers  are  presented  in  Fig¬ 
ures  127-130.  Once  again  it  should  be  remembered  that  Che  field  data  used  to 
compute  the  seasonal  averages  were  very  sparse. 

85.  Figures  131-134  illustrate  a  comparison  of  computed  and  recorded 
stratification  for  the  complete  year  for  the  main  bay  and  tributary  stations. 

A  linear  regression  analysis  shows  the  square  of  the  correlation  coefficient 
to  be  above  0.70  for  both  sets  of  stations.  The  error  in  the  computed  strati¬ 
fication  is  less  than  1  ppt  for  approximately  70  percent  of  the  time. 

86.  One  final  result  is  presented  for  the  1985  simulation.  During  the 


*  Personal  Communication  with  Chesapeake  Bay  Liaison  Office,  US  Environ¬ 
mental  Protection  Agency.  Annapolis,  MD. 
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early  part  of  November,  a  major  storm  moved  over  the  lower  Chesapeake  Bay.  As 
illustrated  in  Figure  135,  this  resulted  in  a  peak  flow  of  nearly  6,000  m^/sec 
on  the  James  River.  Figure  136  reveals  that  the  numerical  model  computed 
virtually  salt-free  surface  waters  for  the  entire  James  River  at  the  flood 
peak.  Although  field  data  were  not  available  at  the  time  of  peak  flow,  it  can 
be  seen  that  the  model  does  an  excellent  job  of  reproducing  salinities  immedi¬ 
ately  before  and  after.  These  results  serve  to  demonstrate  that  the  numerical 
model  does  an  excellent  job  of  reproducing  extreme  events. 

1986  Simulation 


Boundary  conditions 

87.  The  1986  year  was  broken  into  the  following  five  seasons  for  the 
purpose  of  generating  seasonally  averaged  longitudinal  transects  of  salinity: 

Season  1  =»  1  Jan— 16  Feb 
Season  2  =>  17  Feb-3  May 
Season  3  =>  4  May— 2  Aug 
Season  4  =»  3  Aug-9  Nov 
Season  5  =»  10  Nov-31  Dec 

Time  series  results  and  boundary  forcing  data  are  presented  for  a  portion  of 
the  year.  The  ocean  boundary  tide  at  the  Chesapeake  Bay  Tunnel  is  shown  in 
Figure  137.*  Wind  forcing  data  corrected  to  reflect  the  wind  over  open  water 
are  given  in  Figures  138  and  139.  As  in  1984  and  1985,  these  data  are  from 
the  Norfolk  and  BWI  Airports.**  Freshwater  inflows  on  the  James,  Potomac, 
Patuxent,  and  Susquehanna  Rivers  are  given  in  Figure  140  (USGS  1987a  and  b) . 
Surface  heat  exchange  data  for  season  4  are  listed  in  Table  15. 

Results 

88.  Locations  at  which  results  are  presented  are  shown  on  Figure  141. 
Complete  results  are  given  in  Appendix  C.  A  comparison  of  computed  and 
recorded  water  level  at  Solomons,  MD,  is  given  in  Figure  142  with  an  expanded 
plot  of  a  couple  of  tidal  cycles  near  day  201  presented  in  Figure  143.  The 
reason  for  the  "roughness"  in  the  second  computed  ebb  is  unclear  but  is 


*  Personal  Communication  with  National  Ocean  Service,  Asheville,  NC. 

**  Personal  Communication  with  National  Climatic  Data  Center,  Rockville,  MD. 
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probably  related  to  the  rapidly  changing  wind  direction  illustrated  in 
Figure  139. 

89.  As  for  the  previous  2  years,  salinity  comparisons  are  given  for 
four  main  bay  stations  plus  stations  on  the  James,  Potomac,  and  Patuxent 
Rivers.  These  are  shown  in  Figures  144—150.  Results  at  the  shallow  eastern 
shore  station  EE  3.5  are  not  reproduced  quite  as  well  as  in  the  previous 

2  years.  However,  results  at  the  remaining  stations  exhibit  the  same  excel¬ 
lent  reproduction  of  actual  salinity  values  as  well  as  the  observed  stratifi¬ 
cation  patterns.  Difference  measures  are  given  in  Table  12  and  demonstrate 
that  the  verification  for  all  3  years  is  approximately  the  same. 

90.  Temperature  comparisons  at  these  same  stations  are  presented  in 
Figures  151-156.  Generally,  these  comparisons  appear  to  be  about  as  good  as 
previous  results  presented  for  1984  and  1985.  As  reflected  by  the  mean  bias 
error  in  Table  13,  the  computed  temperatures  are  consistently  too  cold.  Note, 
however,  that  temperature  inversions  are  generally  reproduced  well.  Again,  it 
should  be  remembered  that  the  surface  heat  exchange  data  are  computed  using 
meteorological  data  from  the  Patuxent  River  Naval  Station*  and  then  applied 
uniformly  over  the  complete  grid. 

91.  The  longitudinal  plots  of  seasonally  averaged  salinities  for  the 
two  main  bay  transects  and  those  on  the  James  and  Potomac  Rivers  are  given  in 
Figures  157-160.  The  main  bay  transects  are  reproduced  extremely  well.  The 
coarseness  of  the  grid  in  the  tributaries  results  in  the  tributary  transects 
not  comparing  as  well. 

92.  The  final  results  presented  are  the  plots  illustrating  the  ability 
of  the  model  to  reproduce  the  stratification,  as  reflected  by  the  salinity,  of 
the  main  bay  and  tributaries  in  an  overall  average  sense.  These  results  are 
shown  in  Figures  161-164.  The  squares  of  the  correlation  coefficients  from  a 
linear  regression  analysis  of  results  for  both  the  main  bay  and  the  tributar¬ 
ies  are  0.79  and  0.89,  respectively.  The  error  in  the  computed  stratification 
is  generally  less  than  1  ppt  for  70  percent  of  the  time. 


*  Personal  Communication  with  National  Climatic  Data  Center,  Rockville,  MD. 
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PART  VI:  SUMMARY  AND  CONCLUSIONS 


Summary 

93.  A  time-varying  3-D  numerical  hydrodynamic  model  of  Chesapeake  Bay 
has  been  developed,  verified,  and  applied  in  a  production  mode  to  yield  year¬ 
long  flow  fields  to  drive  a  3-D  water  quality  model.  The  hydrodynamic  model, 
CH3D,  makes  computations  on  a  curvilinear  or  boundary-fitted  planform  grid. 
Physical  processes  impacting  baywide  circulation  and  vertical  mixing  that  are 
modeled  include  tides,  wind,  density  effects  (salinity  and  temperature), 
freshwater  inflows,  vertical  turbulence,  and  the  effect  of  the  earth's 
rotation. 

94.  A  successful  verification  of  the  hydrodynamic  model  required  sets 
of  synoptic  data  representing  a  range  of  events.  Three  relatively  extensive 
synoptic  data  sets  were  identified.  Each  of  these  data  sets  was  approximately 
1  month  long.  The  June-July  1980  set  represented  dry  summer  conditions  and 
was  collected  and  provided  by  Dr.  William  Boicourt  of  the  Horn  Point  Environ¬ 
ment  Laboratories  in  Cambridge,  MD.  The  April  1983  and  September  1983  data 
sets  were  collected  by  NOS  and  provided  by  HydroQual,  Inc.  The  April  data 
reflected  a  large  spring  runoff,  whereas  the  September  data  contained  a  wind- 
induced  mixing  event  that  was  aided  by  a  temperature  inversion.  Application 
of  the  model  to  the  June-July  1980  and  April  1983  data  sets  involved  an  early 
version  of  the  model,  whereas  verification  to  the  September  1983  data  set  was 
with  the  final  version  of  the  model.  Verification  of  the  model  consisted  of 
comparing  recorded  and  computed  water  levels,  flow  velocities,  salinities,  and 
temperatures.  These  comparisons  were  made  through  an  analysis  of  time  series 
plots  as  well  as  computed  difference  measures. 

95.  With  the  successful  simulation  of  these  data  sets,  the  3-D  model 
was  then  employed  in  a  production  mode.  Yearlong  simulations  of  1984,  1985, 
and  1986  were  made  to  provide  flow  fields  to  drive  long-term  water  quality 
computations.  These  simulations  also  provided  additional  verification  of  the 
model  through  comparisons  of  computed  tides,  salinities,  and  temperatures  at 
several  interior  stations  with  field  data.  Since  salinity  is  a  conservative 
substance  and  local  concentrations  depend  heavily  upon  the  flow  transport 
process,  if  the  salinity  comparison  is  good,  then  it  is  reasonable  to  assume 
the  flow  is  approximately  correct.  Obviously  diffusion  of  salt  (numerical  or 
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real)  plays  a  part  also,  especially  near  saline  fronts  such  as  those  on  the 
tributaries . 


Conclusions 


96.  Two  basic  conclusions  can  be  drawn.  Based  upon  an  overall  compari¬ 
son  of  model  results  with  available  field  data  from  the  six  simulations 
presented,  it  has  been  demonstrated  that  the  3— D  hydrodynamic  model  behaves 
well  and  is  a  good  representation  of  the  hydrodynamics  of  the  Chesapeake  Bay. 
This  has  been  demonstrated  over  virtually  all  types  of  events  that  occur  in 
the  bay,  e.g. ,  a  200— year  flood  event  on  the  James  River  during  November  1985. 
For  the  three  yearlong  simulations  the  average  relative  error  in  the  computed 
salinity  in  the  main  bay  is  about  10  percent,  whereas  in  the  tributaries  it  is 
15-20  percent.  The  coarseness  of  the  grid  in  the  tributaries  is  obviously  the 
major  reason  for  the  increased  error.  The  average  relative  error  for  the 
temperature  is  slightly  less  than  for  the  salinity.  This  study  has  also 
demonstrated  that  yearlong  3— D  hydrodynamic  computations  to  drive  water  qual¬ 
ity  models  are  feasible.  A  1-year  simulation  on  the  Chesapeake  Bay  grid  using 
a  time-step  of  5  minutes  takes  about  10  hours  on  a  CRAY  Y-MP  computer. 

97.  Other  conclusions  relate  to  the  modeling  approach  taken.  The  major 
problem  encountered  with  the  horizontal  boundary-fitted  transformation  con¬ 
cerned  the  convection  terms  in  the  momentum  equations.  Initially  these  terms 
were  derived  such  that  velocity  squared  terms  were  present.  With  this  form, 
stable  solutions  were  not  possible  on  irregular  grids.  However,  after  deiiv- 
ing  these  terms  such  that  only  gradients  of  velocity  squared  terms  appeared, 
this  problem  was  eliminated.  The  use  of  a  sigma  grid  in  the  vertical  resulted 
in  excess  advection  of  salt  from  the  deep  channels  to  the  shallows.  There¬ 
fore,  to  maintain  stratification  in  the  channels  over  long  periods  of  simula¬ 
tion,  computations  are  now  made  on  the  z-plane. 

98.  The  obvious  benefit  of  generating  solutions  on  boundary-fitted 
grids  is  that  geometric  features  are  modeled  more  accurately  and  economically. 
The  price  for  this  is  the  increased  complexity  of  the  equations.  However, 
since  the  equations  are  differenced  on  a  transformed  regular  grid,  this  poses 
no  particular  problem.  Use  of  the  contravariant  velocity  results  in  a  more 
straightforward  differencing  of  the  equations  since  the  same  type  staggered 
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grid  employed  in  Cartesian  models  can  be  used.  In  addition,  boundary 
conditions  are  more  easily  applied. 
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Table  1 

Svunmarv  of  Conditions  During  Simulation  Periods 


Period 


Jun— Jul 
1980 


Apr  1983 


Sep  1983 


1984 


1985 


1986 


Inflow  Wind  Model 

Conditions  Conditions  Used 


Dry  summer  Normal  Early 

conditions  version 


Spring  Normal  Early 

runoff  version 


Low 


Wind-mixing  Final 
event  version 


High  Wide  range  Final 

version 


Low  Wide  range  Final 

version 


Average  Wide  range  Final 

version 


_ Comments _ 

Data  from  Boicourt.  Over¬ 
land  winds  used  in  computa¬ 
tion  of  surface  heat 
exchange .  Boundary  temp 
and  salt  constant  in  time. 

Data  from  NOS .  Overland 
winds  used  in  computation 
of  surface  heat  exchange . 
Boundary  temp  and  salt 
constant  in  time. 

Data  from  NOS.  Entire  bay 
well  mixed  during  wind 
event.  Over-water  winds 
used  in  surface  heat 
exchange  computations. 

All  boundary  conditions 
time  varying. 

Data  from  NOS  and  USEPA. 
Over-water  winds  used  in 
surface  heat  exchange  com¬ 
putations  .  All  boundary 
conditions  time  varying. 
Little  data  first  6  months. 

Data  from  NOS  and  USEPA. 
Overwater  winds  used  in 
surface  heat  exchange  com¬ 
putations  .  All  boundary 
conditions  time  varying. 
Extreme  inflow  event  in 
November  on  lower  bay. 

Data  from  NOS  and  USEPA. 
Overwater  winds  used  in 
surface  heat  exchange 
computations.  All  boundary 
conditions  time  varying. 


Table  2 

1980  Surface  Heat  Exchange  Data 


Date 

Equilibrium 

Temperature 

°C 

Surface 
Transfer 
Coefficient 
Watts /m^/°C 

6/23/80 

24.7 

58.2 

6/24/80 

24.4 

90.6 

6/25/80 

23.5 

80.2 

6/2f/80 

22.3 

59.8 

6/27/80 

21.7 

31.0 

6/28/80 

24.1 

23.9 

6/29/80 

22.9 

28.6 

6/30/80 

26.9 

49.7 

7/01/80 

27.3 

36.6 

7/02/80 

27.9 

29.8 

7/03/80 

28.0 

32.0 

7/04/80 

26.8 

48.0 

7/05/80 

24.7 

69.1 

7/06/80 

24.1 

47.1 

7/07/80 

22.2 

47.5 

7/08/80 

24.2 

34.4 

7/09/80 

23.4 

47.7 

7/10/80 

23.6 

31.3 

7/11/80 

23.4 

55.7 

7/12/80 

26.7 

43.6 

7/13/80 

25.9 

63.9 

7/14/80 

24.3 

24  8 

7/15/80 

21.3 

38.6 

7/16/80 

22.8 

52.7 

7/17/80 

21.3 

24.6 

7/18/80 

19.5 

33.8 

7/19/80 

21.7 

25.9 

7/20/80 

23.6 

25.3 

7/21/80 

25.4 

28.5 

7/22/80 

27.6 

31.5 

7/23/80 

27.2 

29.2 

7/24/80 

26.6 

40.1 

7/25/80 

25.3 

44.9 

7/26/80 

26.4 

42.8 

7/27/80 

26.3 

45.9 

7/28/80 

27.8 

30.3 

7/2S/80 

25.7 

42.8 

7/30/80 

25.4 

35.2 

7/31/80 

25.5 

43.1 

Table  3 


Date 


4/01/83 

4/02/83 

4/03/83 

4/04/83 

4/05/83 

4/06/83 

4/07/83 

4/08/83 

4/09/83 

4/10/83 

4/11/83 

4/12/83 

4/13/83 

4/14/83 

4/15/83 

4/16/83 

4/17/83 

4/18/83 

4/19/83 

4/20/83 

4/21/83 

4/22/83 

4/23/83 

4/24/83 

4/25/83 

4/26/83 

4/27/83 

4/28/83 

4/29/83 

4/30/83 


April  1983  Surface  Heat  Exchange  Data 


Equilibrium 

Temperature 

°C 

Surface 
Transfer 
Coefficient 
Watts /m^/°C 

13.8 

15. 

7 

9.6 

21. 

7 

14.2 

29. 

4 

16.2 

17. 

4 

14.6 

11. 

8 

13.7 

12. 

4 

14.6 

12. 

8 

16.7 

13. 

6 

12.6 

23. 

0 

13.8 

18. 

7 

12.4 

21. 

9 

18.2 

15. 

3 

15.7 

19. 

2 

13.1 

29, 

,2 

13.6 

33, 

,9 

12.1 

27 

.7 

14.3 

18 

.6 

9.9 

18, 

,  7 

4.6 

26, 

.2 

6.5 

28 

.4 

18.4 

15 

.1 

19.1 

15 

.8 

14.7 

19 

.6 

12.4 

29 

.6 

7.4 

30 

.6 

15.2 

25 

.3 

20.9 

18 

.9 

21.8 

27 

.1 

17.9 

32 

.6 

17.9 

38 

.8 

Table  4 

September  1983  Surface  Heat  Exchange  Data 


Surface 

Equilibrium 

Transfer 

Temperature 

Coefficient 

Date 

°C 

Watts /m^/°C 

9/01/83 

24.2 

41.7 

9/02/83 

23.6 

31.4 

9/03/83 

27.2 

20.0 

9/04/83 

26.0 

45.1 

9/05/83 

28.7 

34.8 

9/06/83 

28.8 

41.4 

9/07/83 

29.4 

27.8 

9/08/83 

23.4 

27.9 

9/09/83 

26.5 

25.1 

9/10/83 

26.9 

29.5 

9/11/83 

25.8 

42.2 

9/12/83 

25.4 

42.9 

9/13/83 

21.3 

43.4 

9/14/83 

17.1 

46.2 

9/15/83 

18.1 

38.1 

9/16/83 

19.4 

29.9 

9/17/83 

22.3 

32.8 

9/18/83 

22.6 

29.6 

9/19/83 

23.4 

49.3 

9/20/83 

23.4 

60.5 

9/21/83 

20.8 

62.9 

9/22/83 

12.9 

46.5 

9/23/83 

15.4 

21.0 

9/24/83 

15.7 

22.3 

9/25/83 

20.6 

13.5 

9/26/83 

18.2 

13.3 

9/27/83 

21.3 

18.7 

9/28/83 

17.5 

40.4 

9/29/83 

14.5 

54.8 

9/30/83 

16.9 

58.3 

9/31/83 

20.8 

22.4 

10/1/83 

19.9 

19.3 

Table  5 

Comparison  of  Model  and  Observed  Water  Surface  at 
Maximum  Flood  and  Minimum  Ebb 


Maximum 

Flood 

Minimum 

Ebb 

Station 

K  -  Eol 

m 

-  tj 
min 

|E„  -  Ej 
m 

1 

min 

Hampton  Roads 

0.026 

12.55 

0.104 

26.42 

Lewisetta 

0.085 

47.50 

0.061 

50.63 

Colonial  Beach 

0.058 

31.25 

0.055 

30.00 

Solomons 

0.113 

48.51 

0.058 

50.00 

Annapolis 

0.085 

15.65 

0.073 

18.51 

Havre  de  Grace 

0.192 

160.85 

0.098 

171.06 

Note;  E„  -  computed  water-surface  elevation  at  maximum  flood  (minimum  ebb) 

Eq  “  observed  water-surface  elevation  at  maximum  flood  (minimum  ebb) 

t„  -  computed  time  to  maximum  flood  (minimum  ebb) 

tg  ^  observed  time  to  maximum  flood  (minimum  ebb) 


Table  6 

Difference  Measures  for  the  September  1983 
Water-Surface  Elevations 


Station 

MBE 

m 

RMSE 

m 

MAE 

m 

Hampton  Roads 

0. 

066 

0. 

095 

0. 

074 

Lewisetta 

0. 

016 

0. 

098 

0. 

077 

Colonial  Beach 

0. 

013 

0. 

,  106 

0. 

.088 

Solomons 

0. 

,067 

0 

.127 

0. 

,  101 

Annapolis 

0 

.026 

0 

.111 

0 

.085 

Havre  de  Grace 

0 

.051 

0 

.377 

0 

.  328 

Table  7 

Comparison  of  Model  and  Observed  Near-Surface  Alonp-Ba 
Current  at  Maximum  Flood  and  Minimum  Ebb 


Station 

Maximum 

Flood 

Minimum 

Ebb 

lv„  -  vj 

cm/sec 

-  t„l 
min 

|v„  -  vj 

cm/sec 

lt„  -  tj 

min 

Wolf  Trap 

10.97 

32.23 

15.81 

35.47 

Mid  Bay 

8.78 

28.61 

10.62 

25.36 

Bay  Bridge 

6.95 

81.84 

7.80 

86.74 

Note;  V„  =  Computed  along-bay  current  near  the  surface  at  maximum  flood 
(minimum  ebb) 

Vq  —  Observed  along-bay  current  near  the  surface  at  maximum  flood 
(minimum  ebb) 

t„  =  Computed  time  to  maximum  flood  (minimum  ebb) 
t„  -  Observed  time  to  maximum  flood  (minimum  ebb) 


Station 

Level 

MBE 

cm/sec 

RMSE 

cm/sec 

MAE 

cm/sec 

Wolf  Trap 

Near  surface 

-2.17 

23.12 

19.18 

Midbay 

Near  surface 

-0.61 

12.03 

9.81 

Near  bottom 

-0.06 

10.73 

8.26 

Bay  Bridge 

Near  surface 

-5.24 

21.43 

18.31 

Near  bottom 

-14.21 

18.97 

16.79 

Table  9 

Difference  Measures  for  the  September  1983  Salinities 


Station 

Level 

MBE 

DPt 

RMSE 

pot 

MAE 

ppt 

ARE 

% 

Wolf  Trap 

Near  surface 

1.28 

2.17 

1.54 

10.22 

Midbay 

Near  surface 

0.26 

0.67 

0.55 

4.57 

Near  bottom 

-0.57 

1.19 

0.87 

5.91 

Bay  Bridge 

Near  surface 

0.14 

0.69 

0.53 

5.15 

Near  bottom 

-0.08 

1.73 

1.03 

8.14 

Table  10 

Difference  Measures  for  the  September  198^  Temperatures 


MBE 

RMSE 

MAE 

ARE 

Station 

Level 

“C 

“C 

°C 

% 

Wolf  Trap 

Near  surface 

-0.60 

0.89 

0.76 

3.60 

Midbay 

Near  surface 

-0.17 

0.48 

0.41 

i  .  94 

Near  bottom 

0.12 

0.38 

0.32 

1.55 

Bay  Bridge 

Near  surface 

0.20 

0.71 

0.52 

2.91 

Near  bottom 

-0.10 

0.21 

0.12 

0.83 

Table  11 

Surface  Heat  Exchange  Data.  Season  5,  1984 


Surface 

Equi''  ibriuiii 

Transfer 

Temperature 

Coefficient 

Date 

°C 

Watts /m^/°C 

11/26/84 

5.3 

13.0 

11/27/84 

5.6 

13.9 

11/28/84 

10.6 

55.0 

11/29/84 

3.3 

32.9 

11/30/84 

4.9 

23.8 

12/01/84 

6 . 4 

27.5 

12/02/84 

4.8 

x5.0 

12/03/84 

8.0 

22.1 

12/04/84 

1.1 

17.6 

12/05/84 

0.5 

25.2 

12/06/84 

1.1 

36.6 

12/07/84 

-4.4 

35.0 

12/08/84 

0.3 

20.2 

12/09/84 

3.8 

10.0 

12/10/84 

3.8 

10.3 

12/11/84 

5.3 

19.0 

12/12/84 

4.1 

17.1 

12/13/84 

7.8 

14  1 

12/14/84 

8.6 

14.4 

12/15/84 

8.5 

15.4 

12/16/84 

6.7 

17.1 

12/17/84 

10.5 

20.1 

12/18/84 

13.0 

16.4 

12/19/84 

9.1 

19.1 

12/20/84 

6.8 

10.8 

12/21/84 

6.1 

22.0 

12/22/84 

11.4 

53.3 

:'>/23/84 

2.3 

20.6 

12/24/84 

4.0 

21.5 

12/25/84 

3.4 

33.5 

12/26/84 

1.0 

9.3 

12/27/84 

3.5 

11.9 

12/28/84 

11.4 

30.6 

12/29/84 

13.6 

51.5 

12/30/84 

9.9 

30.2 

12/31/84 

9.5 

31.1 

Table  12 

Salinity  Statistics  for  Yearlong  Simulations 


Main 

Bay 

Tributaries 

MBE 

MAE 

RMSE 

ARE 

MBE 

MAE 

RMSE 

ARE 

Year 

ppt 

ppt 

ppt 

% 

ppt 

ppt 

% 

r’::ar  Surface 

1984 

-0.05 

1.21 

1./.5 

10.13 

-0.20 

1.54 

1.85 

17.25 

1985 

0.90 

1.49 

1.93 

11.43 

1.09 

1.66 

2.04 

16.32 

1986 

0.25 

1.52 

1.87 

12.29 

0.32 

1.60 

2.04 

15.02 

Near  Bottom 

1984 

-0.64 

1.72 

2.12 

10.44 

-0.09 

2.36 

2.83 

19.94 

1985 

1.12 

1.78 

2.34 

12.11 

1.27 

2.31 

2.99 

18.30 

1986 

-0.16 

1.35 

1.68 

8.74 

0.24 

2.23 

2.65 

13.40 

Table  13 

Temperature  Statistics  for  Yearlone  Simulations 


Main 

Bay 

Tributaries 

MBE 

MAE 

RMSE 

ARE 

MBE 

MAE 

RMSE 

ARE 

Year 

°C 

°C 

'’C 

% 

°C 

°C 

°C 

% 

Near  Surface 

1984 

-1.48 

1.56 

1.79 

8.45 

-1.90 

2.07 

2.20 

11.06 

1985 

-1.73 

1.97 

2.25 

12.76 

-2.03 

2.06 

2.28 

12.11 

1986 

-1.61 

1.72 

2.04 

11.62 

-1.97 

2.09 

2.40 

14.92 

Near  Bottom 

1984 

-0.42 

1.14 

1.35 

6.79 

-1.28 

1.70 

1.92 

9.70 

1985 

-1.01 

1.45 

1.70 

9.99 

-1.53 

1.75 

2.00 

10.94 

1986 

-1.03 

1.34 

1.63 

9.60 

-1.61 

1.91 

2.32 

14.74 

Table  14 

Surface  Heat  Exchange  Data.  Season  4.  1985 


Date 

Equilibrium 

Temperature 

°C 

Surface 
Transfer 
Coefficient 
Watts /m^/°C 

7/16/85 

27.1 

35.8 

7/17/85 

24.0 

40.2 

7/18/85 

29.9 

21 . 3 

7/19/85 

25.5 

38.4 

7/20/85 

27.6 

32 . 6 

7/21/85 

27.4 

29.0 

7/22/85 

25.1 

43 . 1 

7/23/85 

22.2 

34.3 

7/24/85 

24.8 

31.4 

7/25/85 

24.6 

65 . 0 

7/26/85 

24.2 

88.8 

7/27/85 

25.0 

26 . 2 

7/28/85 

27.5 

24.9 

7/29/85 

26.4 

17.8 

7/30/85 

26.7 

38.2 

7/31/85 

25.5 

53.2 

8/01/85 

21.9 

40 . 1 

8/02/85 

22.5 

33.7 

8/03/85 

27.9 

18.9 

8/04/85 

24.9 

27 . 5 

8/05/85 

23.2 

32.8 

8/06/85 

23.1 

28 . 7 

8/07/85 

23.2 

56 . 0 

8/08/85 

25.3 

30.1 

8/09/85 

27.8 

28 . 3 

8/10/85 

26.1 

35.7 

8/11/85 

30.7 

20.4 

8/12/85 

24.6 

33 . 1 

8/13/85 

25.3 

44 . 2 

8/14/85 

26.3 

44 . 2 

3/15/85 

24.9 

60.5 

8/16/85 

27.1 

33 . 5 

8/17/85 

22.0 

35 . 6 

8/18/85 

21.2 

63 . 1 

8/19/85 

24.3 

49 . 8 

8/20/85 

27.5 

20.9 

8/21/85 

21.9 

35 . 5 

8/22/85 

23.0 

22 . 1 

8/23/85 

26.2 

17.4 

8/24/85 

19.5 

54 . 5 

8/25/85 

21.8 

50.8 

8/26/85 

24.6 

44 . 5 

8/27/85 

23.6 

33 . 2 

(Continued) 


Table  14  (Concluded) 


Date 

Equilibrium 

Temperature 

°C 

Surface 
Transfer 
Coefficient 
Watts /m^/°C 

8/28/85 

23.6 

37.5 

8/29/85 

23.0 

53.3 

8/30/85 

22. h 

74.2 

8/31/85 

20.1 

55.3 

9/01/85 

20.5 

56.1 

9/02/85 

22.9 

57.3 

9/03/85 

26.0 

36.2 

9/04/85 

26.1 

41.5 

9/05/85 

28.1 

34.2 

9/06/85 

27.0 

44.7 

9/07/85 

30.7 

21.1 

9/08/85 

28.9 

20.3 

9/09/85 

27.0 

30.4 

9/10/85 

28.2 

27.6 

9/11/85 

20.6 

35.0 

9/12/85 

17.4 

28.5 

9/13/85 

13.0 

43.3 

9/14/85 

16.3 

27.3 

9/15/85 

18.7 

19.3 

9/16/85 

20.4 

19.1 

9/17/90 

20.8 

21.7 

9/18/90 

22.9 

15.0 

9/19/90 

22.5 

20.5 

Table  15 

Surface  Heat  Exchange  Data.  Season  A.  1986 


Date 

Equilibrium 

Temperature 

°C 

Surface 
Transfer 
Coefficient 
Watts /m^/°C 

8/03/86 

24.7 

33.0 

8/04/86 

26.3 

30.2 

8/05/86 

26.5 

34.1 

8/06/86 

23.9 

53.6 

8/07/86 

26.0 

36.9 

8/08/86 

24.3 

50.7 

8/09/86 

25.8 

39.2 

8/10/86 

25.2 

56.1 

8/11/86 

22.9 

52.3 

8/12/86 

20.5 

34.3 

8/13/86 

20.2 

50.3 

8/14/86 

?4.2 

42.4 

8/15/86 

26.2 

42.4 

8/16/86 

26.6 

55.6 

8/17/86 

25.5 

60.6 

8/18/86 

23.7 

72.8 

8/19/86 

24.5 

32.7 

8/20/86 

21.4 

83.9 

8/21/86 

22.7 

56.0 

8/22/86 

23.8 

32.0 

8/23/86 

22.6 

52.7 

8/24/86 

20.1 

62.2 

8/25/86 

21.7 

26.3 

8/26/86 

23.3 

50.4 

8/27/86 

22.9 

86.7 

8/28/86 

16.3 

55.9 

8/29/86 

16.9 

30.5 

8/30/86 

18.2 

31.2 

8/31/86 

20.9 

64.3 

9/01/86 

20.9 

51.9 

9/02/86 

19.7 

33.6 

9/03/86 

20.9 

57.3 

9/04/86 

21.3 

73.3 

9/05/86 

21.4 

34.8 

9/06/86 

24.3 

39.9 

9/07/86 

24.5 

20.1 

9/08/86 

17.3 

29.8 

9/09/86 

19.2 

23.1 

9/10/86 

19.5 

47.4 

9/11/86 

22.0 

50.4 

(Continued) 


(Sheet  1  of  3) 


Table  15  (Continued) 


Date 

Equilibrium 

Temperature 

°C 

Surface 
Transfer 
Coefficient 
Watts /m^/°C 

9/12/86 

22.2 

31.5 

9/13/86 

19.5 

29.8 

9/14/86 

20.5 

39.7 

9/15/86 

20.1 

34.2 

9/16/86 

15.8 

19.3 

9/17/86 

13.2 

17.9 

9/18/86 

15.1 

35.1 

9/19/86 

19.5 

30.0 

9/20/86 

19.9 

21.7 

9/21/86 

21.6 

23.5 

9/22/86 

19.3 

22.2 

9/23/86 

21.8 

26.6 

9/24/86 

22.9 

27.2 

9/25/86 

24.8 

40.2 

9/26/G5 

25.6 

39.6 

9/27/86 

21.6 

25.4 

9/28/86 

19.2 

34.9 

9/29/86 

22.4 

33.2 

9/30/86 

23.6 

23.1 

10/01/86 

23.9 

27.1 

10/02/86 

23.7 

39.1 

10/03/86 

22.2 

47.0 

10/04/86 

23.0 

40.4 

10/05/86 

20.6 

50.7 

10/06/86 

12.7 

28.6 

10/07/86 

14.0 

20.9 

10/08/86 

17.7 

39.5 

10/09/86 

21.0 

26.1 

10/10/86 

10.1 

45.1 

10/11/86 

12.2 

105.3 

10/12/86 

15.4 

97.2 

10/13/86 

17.9 

32.4 

10/14/86 

18.4 

29.6 

10/15/86 

8.8 

44.7 

10/16/86 

12.1 

61.5 

10/17/86 

11.4 

39.8 

10/18/86 

10.4 

50. 7 

10/19/86 

10.3 

53.7 

10/20/86 

11.7 

48.5 

10/21/86 

12.5 

47.2 

10/22/86 

15.1 

48.7 

10/23/86 

16.7 

82.3 

10/24/86 

14.3 

55.8 

(Continued) 
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Table  15  (Concluded) 


Surface 


Equilibrium 

Transfer 

Temperature 

Coefficient 

Date 

°C 

Watts /m^/°C 

10/25/86 

11.7 

32.5 

10/26/86 

16.1 

33.4 

10/27/86 

16.8 

48.5 

10/28/86 

13.7 

34.2 

10/29/86 

12.9 

36.1 

10/30/86 

13.5 

33.7 

10/31/86 

9.3 

37.6 

11/01/86 

14.6 

21.9 

11/02/86 

14.7 

28.6 

11/03/86 

8.1 

26.1 

11/04/86 

10.9 

25.9 

11/05/86 

9.3 

31.7 

11/06/86 

9.4 

30.4 

11/07/86 

10.8 

21.8 

11/08/86 

16.2 

32.6 

11/09/86 

16.7 

65.4 

(Sheet  3  of  3) 


Figure  1 .  Chesapeake  Bay 


Figure  2.  Planform  boundary-fitted  grid  of  Chesapeake  Bay 
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Figure  4.  Ocean  boundary  tide  at  the  Chesapeake  Bay  Tunnel  during 

June-July  1980 
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Figure  5.  Freshwater  inflows  during  June-July  1980. 
Day  0  is  23  June.  (Sheet  1  of  3) 
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Figure  5.  (Sheet  2  of  3) 
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g.  Susquehanna  River 
Figure  5.  (Sheet  3  of  3) 
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Figure  8.  Comparison  of  computed  and  recorded  tide  at  Hampton 
Roads,  VA,  during  June-July  1980.  Day  0  is  23  June. 
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Figure  9.  Comparison  of  computed  and  recorded  tide  at  Solomons, 
•ID,  during  June-July  1980.  Day  0  is  23  June. 
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Figure  10.  Comparison  of  computed  and  recorded  tide  at  Annapolis, 
MD,  during  June-July  1980.  Day  0  is  23  June. 
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Figure  11.  Comparison  of  computed  and  recorded  tide  at  Havre  dc 
Grace,  MD,  during  June-July  1980.  Day  0  is  23  June. 


b.  Depth  -  27.5  ft 


Figure  13.  Comparison  of  computed  and  recorded 
velocity  at  station  WT4  during  June-July  1980. 
Day  0  is  23  June. 
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Figure  14.  Comparison  of  computed  and  recorded  salinity  at 
station  WT4  during  June-July  1980.  Day  0  is  23  June. 
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c.  Depth  =  57.5  ft 


Figure  15.  Comparison  of  computed  and  recorded 
temperature  at  station  WT4  during  June-July  1980. 
Day  0  is  23  June. 


a.  Depth  -  7.5  ft 


b.  Depth  —  22.5  ft 

Figure  16 .  Comparison  of  computed  and  recorded 
velocity  at  station  SP2  during  June-July  1980. 
Day  0  is  23  June. 


Figure  17.  Comparison  of  computed  and  recorded 
salinity  at  station  SP2  at  a  depth  of  12.5  ft 
during  June-July  1980.  Day  0  is  23  June. 
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Figure  19.  Comparison  of  computed  and  recorded 
velocity  at  station  BBl  at  a  depth  of  12.5  ft 
during  June-July  1980.  Day  0  is  23  June. 
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Figure  20.  Comparison  of  computed  and  recorded  salinity  at 
station  BBl  during  June-July  1980.  Day  0  is  23  June. 


TEMPERATURE.  °C  TEMPERATURE. 


o 


o  i 

C>4  ' 


I  LEGEND 
COMPUTED 
'  OBSERVHT 


Of . . . . . . .  . . . 

5  10  15  20  25  30  35  40 

a.  Depth  ^  12.5  ft 


o 

n 


CN 


o 

CN 


LEGEND 

COMPUTED 

OBSERVED 


o  . 


5 


10  15  20  25  30  35 

b.  Depth  =  52.5  ft 


40 


Figure  21,  Comparison  of  computed  and  recorded  temperature  at 
station  BBl  during  June-July  1980.  Day  0  is  23  June. 
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Figure  22.  Location  of  data  stations  for  April  1983  data  set 
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Figure  23.  Wind  data  at  Patuxent  River  Naval  Air  Station 
during  April  1983.  Day  0  is  1  April. 


Figure  24.  Wind  data  << rfolk  International  Airport  during 
April  1^63  Day  0  is  1  April. 
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Figure  23.  Freshwater  inflows  during  April 
1983.  Dav  0  is  1  April.  (Sheet  1  of  3) 
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Figure  25.  (Sheet  2  of  3) 
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Figure  26.  Ocean  boundary  tide  at  Chesapeake  Bay  Tunnel  during 

April  1983 
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Figure  27.  Comparison  of  computed  and  recorded  tides  at  Hampton 

Roads,  VA,  during  April  1983 
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Figure  28. 


Comparison  of  computed  and  recorded  tides  at  Colonial 
Beach,  VA,  during  April  1983 
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Figure  29.  Comparison  of  computed  and  recorded  tides  at  Solomons,  MLV 

during  April  1983 
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Figure  30.  Comparison  of  computed  and  recorded  tides  at  Annapolis 

MD,  during  April  1983 
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Figure  31.  Comparison  of  computed  and  recorded  tides  at 
Havre  de  Grace  during  April  1983 
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Figure  32.  Comparison  of  computed  and  recorded  velocity  at 
station  BE  during  April  1983 
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Figure  33.  Comparison  of  computed  and  recorded  velocity  at 
station  WT  during  April  1983 
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Figure  34.  Comparison  of  computed  and  recorded  velocity  at 
station  MB  during  April  1983 
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b.  Depth  “  42.5  ft 


Figure  35.  Comparison  of  computed  and  recorded  salinity  at 
station  MB  during  April  1983 
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Figure  36.  Comparison  of  computed  and  recorded  temperature  at 
station  MB  during  April  1983 
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b.  Depth  -  37.5  ft 
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c.  Depth  -  62.5  ft 


Figure  37.  Comparison  of  computed  and 
recorded  velocity  at  station  BB  during 
April  1983 
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a.  Depth  =  22.5  ft 
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b.  Depth  =  37.5  ft 
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Figure  38.  Comparison  of  computed  and 
recorded  salinity  at  station  BB  during 
April  1983 
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Figure  39.  Comparison  of  computed  and 
recorded  temperature  at  station  BB 
during  April  1983 


Baltimore 
Washington 
Airport 


DISCHARGE,  CMS  DISCHARGE,  CMS  DISCHARGE,  CMS 


Figvire  h''  Freshwater  inflows  during  September  1983.  Day  0  is 

1  September.  (Sheet  1  of  3) 
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Figure  44.  (Sheet  2  of  3) 
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Figure  44.  (Sheet  3  of  3) 
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Figure  45.  Wind  data  at  the  Norfolk  International  Airport  during 
September  1983.  Day  0  is  1  September. 
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Figure  46.  Wind  data  at  the  Patuxent  River  Naval  Air  Station  during 
September  1983.  Day  0  is  1  September. 
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Figure  48.  Comparison  of  computed  and  recorded  tide  at  Lewisetta 

VA,  during  September  1983 
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Figure  49.  Comparison  of  computed  and  recorded  tide  at  Colonial 
Beach,  VA,  during  September  1983 
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Figure  50.  Comparison  of  computed  and  recorded  tide  at  Solomons, 

MD,  during  September  1983 
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Figure  51.  Comparison  of  computed  and  recorded  tide  at  Annapolis, 

MD,  during  September  1983 
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Figure  52.  Comparison  of  computed  and  recorded  tide  at  Havre 
de  Grace,  MD,  during  September  1983 
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Figure  53.  Filtered  ocean  boundary  tide  during  September  1983 
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Figure  54. 
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Comparison  of  filtered  computed  and  recorded  tide  at 
Hampton  Roads,  VA,  during  September  1983 
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Figure  55.  Comparison  of  filtered  computed  and  recorded  tide  at 
Colonial  Beach,  VA,  during  September  1983 


Ml 


(O 

QC 

Ui 


OC  o  '  i 

iu  '  '  LEGEND 

St  COMPUTED 

i ^  OBSERVED 

•  I  i 


Ml  , 

7  14  21  28 

Figure  56.  Comparison  of  filtered  computed  and  recorded  tide  at 
Lewisetta,  VA,  during  September  1983 
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Figure  57 .  Comparison  of  filtered  computed  and  recorded  tide  at 
Solomons,  MD,  during  September  1983 
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Figure  58.  Comparison  of  filtered  computed  and  recorded  tide  at 
Annapolis,  MD,  during  September  1983 


Figure  60.  Comparison  of  computed  and  recorded  velocity  at 
station  WT  during  September  1983 


200  -100  0  100  200 


7  14  21  28 


b .  Depth  -42.5  ft 

Figure  61.  Comparison  of  computed  and  recorded  velocity  at 
station  MB  during  September  1983 
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Figure  63.  Comparison  of  computed  and  recorded  salinity  at 
station  WT  during  September  1983 
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Figure  64.  Comparison  of  computed  and  recorded  salinity  at 
station  MB  during  September  1983 
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Figure  65.  Comparison  of  computed  and  recorded  salinity  at 
station  BB  during  September  1983 
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Figure  68 .  Comparison  of  computed  and  recorded  temperature  at 
station  BB  during  September  1983 
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Figure  69.  Comparison  of  computed  and  recorded  salinity  at  a  depth  of 
42.5  ft  at  station  MB  without  the  three-point  smoothing 
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Figure  73.  Eudy  diffusiv 
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Figure  73.  (Sheet  3  of  4) 
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Figure  76.  Wind  data  at  the  Norfolk  International  Airport  during  1984. 
Day  0  corresponds  to  day  262  of  1984. 
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Figure  79.  Data  stations  for  1984  data  set 
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Figure  80.  Comparison  of  computed  and  recorded  tide  at 
Annapolis,  MD,  during  1984 
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Figure  81.  Comparison  of  computed  and  recorded  tide  at 
Annapolis,  MD,  near  day  301  of  1984 
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re  82.  Comparison  of  computed  and  recorded  salinity  at 
station  EE  3.5  during  1984 


Figure  83.  Comparison  of  computed  and  recorded  salinity  at 
station  CB  6.3  during  1984 
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84.  Comparison  of  computed  and  recorded  salinity  at 
station  CB  5.1  during  season  4  of  1984 
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Figure  87 .  Comparison  of  computed  and  recorded  salinity  at 
station  LE  5.5  during  1984 
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e  83.  Comparison  of  computed  and  recorded  salinity  at 
station  LE  2.2  during  1984 
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Figure  89.  Comparison  of  computed  and  recorded  salinity  at 
station  LE  1.1  during  1984 


Figure  90.  Comparison  of  computed  and  recorded  temperature 
at  station  EE  3.5  during  1984 
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Figure  91.  Comparison  of  computed  and  recorded  temperature 
at  station  CB  6.3  during  1984 
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Figure  93.  Comparison  of  computed  and  recorded  temperature  at 

station  CB  3.3W  during  1984 


DAY  IN  1984 


Figure  94.  Comparison  of  computed  and  recorded  temperature 
at  station  LE  5.5  during  1984 
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Figure  95.  Comparison  of  computed  and  recorded  temperature 
at  station  LE  2.2  during  1984 


Figure  96.  Comparison  of  computed  and  recorded  temperature 
at  station  LE  1.1  during  1984 
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Figure  97.  Comparison  of  seasonally  averaged  salinities  along 

the  main  bay  during  1984 
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98.  Comparison  of  seasonally  averaged  salinities  along 
the  eastern  main  bay  transect  during  1984 
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Figure  99.  Comparison  of  seasonally  averaged  salinities  along 

the  James  River  during  1984 
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Figure  100.  Comparison  of  seasonally  averaged  salinities  along 
the  Potomac  River  during  1984 
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OBSERVED  STRATIFICATION 

Figure  101.  Comparison  of  computed  and 
recorded  stratification  for  all  main 
bay  stations  during  1984 
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Figure  102 .  Frequency  of  occurrence  of  the 
error  in  the  computed  stratification  for  all 
main  bay  stations  during  1984 
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Figure  103.  Comparison  of  computed  and 
recorded  stratification  for  all  tributary 
stations  during  1984 
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Figure  104.  Frequency  of  occurrence  of  the 
error  in  the  computed  stratification  for  all 
tributary  stations  during  1984 
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Figure  105.  Ocean  boundary  tide  at  the  Chesapeake  Bay  Tunnel  during  1985 


Figure  106.  Wind  data  at  the  Norfolk  International  Airport  during  1985. 

Day  0  corresponds  to  Day  195  of  1985. 
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Figure  107.  Wind  data  at  the  BWI  Airport  during  1985.  Day  0  corresponds  to 

Day  195  of  1985. 
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figure  110.  Comparison  of  computed  and  recorded  tide  at 
Colonial  Beach,  VA,  during  1985 
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Figure  111.  Comparison  of  computed  and  recorded  tide 
at  Colonial  Beach,  VA,  near  day  201  during  1985 
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Figure  112.  Comparison  of  computed  and  recorded  salinity 
at  station  EE  3.5  during  1985 
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Figure  113.  Comparison  of  computed  and  recorded  salinity 
at  station  CB  6.3  during  1985 
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Figure  114.  Comparison  of  computed  and  recorded  salinity 
at  station  CB  5.1  during  1985 
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Figure  115.  Comparison  of  computed  and  recorded  salinity  at 
station  CB  3.3W  during  1985 
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Figure  116.  Comparison  of  computed  and  recorded  salinity 
at  station  LE  5.5  during  1985 
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Figure  117.  Comparison  of  computed  and  recorded  salinity 
at  station  LE  2.2  during  1985 
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Figure  118.  Comparison  of  computed  and  recorded  salinity 
at  station  LE  1.1  during  1985 
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Figure  119.  Comparison  of  computed  and  recorded  temperature 
at  station  EE  3.5  during  1985 
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Figure  129.  Comparison  of  computed  and  recorded  temperature 
at  station  CB  6.3  during  1985 
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Figure  121.  Comparison  of  computed  and  recorded  temperature 
at  station  CB  5.1  during  1985 
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Comparison  of  computed  and  recorded  temperature 
at  station  CB  3.3W  during  1985 
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Figure  123.  Comparison  of  computed  and  recorded  temperature 
at  station  LE  5.5  during  1985 
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;ure  124.  Comparison  of  computed  and  recorded  temperature 
at  station  LE  2.2  during  1985 
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Figure  125.  Comparison  of  computed  and  recorded  temperature 
at  station  LE  1.1  during  1985 
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Figure  127.  Comparison  of  seasonally  averaged  salinities  along 

the  main  bay  during  1985 
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Figure  128.  Comparison  of  seasonally  averaged  salinities  along 
the  eastern  main  bay  transect  during  1985 
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OBSERVED  STRATIFICATION 

Figure  131.  Comparison  of  computed  and 
recorded  stratification  for  all  main  bay 
stations  during  1985 
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Figure  132.  Frequency  of  occurrence  of  the 
error  in  the  computed  stratification  for  all 
main  bay  stations  during  1985 
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Figure  133.  Comparison  of  computed  and 
recorded  stratification  for  all  tributary 
stations  during  1985 


IPREDICTED-OBSERVEDI 


Figure  134.  Frequency  of  occurrence  of  the 
error  in  the  computer  stratification  for  all 
tributary  stations  during  1985 
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Figure  135.  Freshwater  inflow  on  James  River  during  1985 

storm  event 
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Figure  136.  Comparison  of  computed  and  recorded 
salinity  at  station  LE  5.5  during  1985  storm 

event 
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Figure  139.  Wind  data  at  the  BWI  Airport  during  1986.  Day  0  corresponds  to 

Day  195  of  1986. 
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Figure  140.  Freshwater  inflows  during  1986  (Continued) 
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Figure  141.  Data  stations  for  the  1986  data  set 
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Figure  144.  Comparison  of  computed  and  recorded  salinity  at 
station  EE  3.5  during  1986 
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Comparison  of  computed  and  recorded  salinity  at 
station  CB  6,3  during  1986 
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Figure  146.  Comparison  of  computed  and  recorded  salinity  at 
station  CB  5.1  during  1986 
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147.  Comparison  of  computed  and  recorded  salinity  at 
station  CB  3 . 3W  during  1986 
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Figure  148.  Comparison  of  computed  and  recorded  salinity  at 
station  LE  5.5  during  1986 
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Figure  149.  Comparison  of  computed  and  recorded  salinity  at 
station  LE  2.2  during  1986 
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Figure  150.  Comparison  of  computed  and  recorded  salinity  at 
station  LE  1.1  during  1986 
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;ure  153.  Compar  on  of  computed  and  recorded  temperature 
at  s  .ation  CB  3 . 3W  during  1986 
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Figure  154.  Comparison  of  computed  and  recorded  temperature 
at  station  LE  5.5  during  1986 
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'ure  155.  Comparison  of  computed  and  recorded  temperature 
at  station  LE  2.2  during  1986 
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Figure  156.  Comparison  of  computed  and  recorded  temperature 
at  station  LE  1.1  during  1986 
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Figure  157. 
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Comparison  of  seasonally  averaged  salinities 
along  the  main  bay  during  1986 
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Figure  158.  Comparison  of  seasonally  averaged  salinities 
along  the  eastern  main  bay  transect  during  1986 
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Figure  159.  Comparison  of  seasonally  averaged  salinities 


along  the  James  River  during  1986 
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Figure  160.  Comparison  of  seasonally  averaged  salinities 
along  the  Potomac  River  during  1986 
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OBSERVED  STRATIFICATION 

Figure  163.  Comparison  of  computed  and 
recorded  stratification  for  all  tributary 
stations  during  1986 
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Figure  164.  Frequency  of  occurrence  of  the 
error  in  the  computed  stratification  for  all 
tributary  stations  during  1986 


APPENDIX  D;  NOTATION 


K.K 

ARE 

C 


Cd 

D 

f 


g 

h 

H 

k 

K 

MAE 

MBE 

Mi 

N 

Oi 

R1 

RMSE 

S 

t 

T 

T, 

u,v 
u'  ,v' ,  w' 
u,  V 
U,V 

Ui.vi 

w 

W 


X ,  y  ,  z 
Zo 

r 

A 


I  ,»7 


Vertical  turbulent  eddy  coefficients 
Average  relative  error 

Surface  drag  coefficient;  conductivity 

Bottom  friction  coefficient 

Eddy  diffusivity  or  viscosity 

Coriolis  parameter 

Acceleration  due  to  gravity 

Layer  thickness 

Total  water  depth 

von  Karman  constant 

Surface  heat  exchange  coefficient 

Mean  absolute  error 

Mean  bias  error 

Model  results 

Brunt-Vaisala  frequency;  number  of  data  points 

Observed  data 

Richardson  number 

Root  mean  square  error 

Salinity 

Time 

Temperature 

Equilibrium  temperature 

Horizontal  components  of  Cartesian  velocity 
Turbulent  velocity  fluctuations 
Components  of  contravar iant  velocity 
Vertically  integrated  contravar iant  unit  flows 

Horizontal  contravar iant  velocity  components  next  to  the  bottom 
Vertical  component  of  velocity 
Wind  speed 

Cartesian  coordinates 
Bottom  roughness  height 
Water-surface  elevation 
Macroscale  of  turbulence 
Boundary-fitted  coordinates 


D1 


Water  density 


Components  of  bottom  shear  stress 
Components  of  surface  shear  stress 


